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1 Introduction 

Photon production has long been considered a key "hard probe" for studying the formation 
and evolution of the quark-gluon plasma in heavy ion collisions. A chief advantage is that 
the coupling of the plasma to photons is weak, which means that the re-absorption rate of 
photons is expected to be negligible. Once formed, a photon will escape to the detector, 
carrying direct information about its formation process unmodified by hadronization or 
other late time physics. 

Experimentally, there are now detailed data on real photon production at RHIC [1-3] 
and the LHC [4-7]. Photons arising from meson decays following hadronization are sub- 
tracted from the data experimentally, and the remaining sample of direct photons arises 
from several (hopefully distinct) processes. There are "prompt" photons produced in the 
scattering of the partons from the colliding nuclei. The production rate here should be 
calculable using perturbative QCD [8]. There are also photons associated with the frag- 
mentation of jets and with jet-medium interactions [9, 10], and photons produced by the 
interaction of excitations of the nearly thermal Quark-Gluon Plasma, which appears to be 
produced in the collision. The thermal and jet-medium photons are the most interesting 
(to us) because they represent a signal specifically of the plasma and its evolution. 

On the theoretical side, the calculation of the photon production rate from the quark- 
gluon plasma has mostly been carried out within the context of the perturbative or weak- 
coupling expansion. The photon emission rate from the plasma was computed to leading 
order in the logarithm of the strong coupling in 1991, when Kapusta et al and Baier et al 
computed the rate of Compton and pair annihilation processes [11, 12]. It was later pointed 
out that this calculation is not complete at leading order, as bremsstrahlung processes arise 
at the same power of as [13] . The complete treatment of these processes was completed in 
2001, when Arnold, Moore and Yaff'e performed a leading-order calculation of the photon 
production rate from an equilibrium plasma [14, 15]. All these calculations are for the 
production of photons from a thermal medium, but they are rather easily adapted to 
include jet-medium interaction photons as well [10]. 

A leading order calculation begs many questions. By itself it gives no information on 
its own reliability; we do not know how quickly the perturbative expansion will converge, 
or what will be the sign of the next correction. There is some concern that a leading order 
calculation will not be very reliable for the photon production rate. First of all, the coupling 
is rather large at the modest temperatures achieved in heavy ion collisions. Second, the 
convergence of other perturbative expansions in the context of finite temperature QCD 
is not very comforting. For instance, the pressure of the Quark-Gluon Plasma has been 
computed to very high order in the perturbative expansion [16-19], with the result that 
the convergence of the series is rather poor. The pressure is a thermodynamical quantity, 
as in fact are almost all the quantities computed beyond leading order in the coupling. 
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Among dynamical transport quantities, such as the photon rate and the rather closely 
related shear viscosity, heavy quark energy loss, and heavy quark diffusion rate, only one 
quantity is known beyond leading order - the heavy quark diffusion rate [20]. In this case 
the next-to-leading order corrections prove to be very large. But this case also may not 
be very representative, since it involves rather different physics than photon production 
or shear viscosity. Heavy quark diffusion involves particles which are nearly at rest and 
interact only via spacelike longitudinal gluons. But the other transport coefficients involve 
species moving at almost the speed of light, exchanging transverse and longitudinal gluons 
at finite frequencies. They also involve light quarks and their hard thermal loops in a much 
more direct way than the heavy quark calculation did. 

We therefore think it would be extremely useful to compute the photon production 
rate at next-to-leading order in the coupling. As we have emphasized, this may be of phe- 
nomenological interest. And it is most definitely of theoretical interest, since it extends our 
understanding of the convergence properties of the perturbative expansion for dynamical 
quantities in the Quark-Gluon Plasma. It also allows us to develop the theoretical tools 
and understanding which will be needed for other quantities, such as the shear viscosity, 
in the context of a very cleanly defined calculation. In the remainder of this paper we will 
present precisely this next-to-leading order calculation of the photon production rate from 
an equilibrium quark-gluon plasma. 

The reason that there are 0{g) NLO corrections is because of the complex self- 
interactions of soft highly occupied gauge fields. As such, all of the NLO corrections 
arise from the interaction with soft gluons. Since we are concerned with hard photons with 
k > T, the photons must be produced by hard quarks moving at essentially the speed of 
light. But such quarks only "see" the soft fields in an eikonalized way, feeling soft-sector 
correlations at lightlike separated points. And bosonic soft correlators at spacelike or light- 
like separated points can be determined from correlators of the Euclidean theory; in fact, 
at leading and next-to-leading order, they are the correlators of the 3-D Euclidean theory, 
EQCD [21]. Better still, the hard particles are only sensitive, at NLO, to two effects from 
the soft physics: a shift in the dispersion relation, and transverse momentum exchange. 
And both of these properties are already known at NLO [21, 22]. 

The goal of our paper is to derive and explain these facts within the context of the 
calculation of the NLO photon production rate. We begin with an overview of the calcu- 
lation in Section 2. We start by reviewing the leading-order calculation, which arises from 
two distinct kinematic regions, one running from "hard" to "soft" fermionic momenta and 
one involving "collinear" fermionic momenta. Next we show that the first corrections to 
the photon production rate arise at 0{g), not at 0{g'^), and that these corrections arise 
in several kinematic regions; the collinear region, the "semi-collinear" region which lies 
between the hard and collinear regions, and the "soft" infrared region. We handle the 
collinear region in Section 3, the soft region in Section 4, and the "semi-collinear" region 
in Section 5. We present our results in Section 6. 

Because the hard modes experience the soft modes at lightlike separations, we can use 
analyticity and Euclidean methods very successfully in the calculation. This was impos- 
sible in the heavy quark diffusion calculation [20], and means that the photon production 
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calculation is actually simpler than the heavy quark diffusion calculation. This technical 
development makes it more likely that other transport coefficients, such as shear viscosity, 
can be computed beyond leading order with a reasonable amount of effort. 

On the phenomenological side, we find that there are rather large corrections to the 
photon production rate from the different kinematic regions, but that they are of both sign 
and surprisingly similar magnitude. In practice the partial NLO corrections nearly cancel 
in the full result for the phenomenologically interesting energy range of several times the 
temperature. As far as we can tell this is an accident. However, the individual (canceling) 
contributions are also not as large, relative to the leading-order result, as in the heavy 
quark case. So it appears plausible that NLO corrections are in general not as severe as 
was feared. 

2 Overview of the calculation 

The photon production rate is given at leading order in a by 



where ea(/c) are a basis of transverse polarization vectors and VF^(i^) is the backward 
Wightman correlator of the electromagnetic current 



Here and throughout the paper capital letters stand for four-vectors, lowercase italic letters 
for the modulus of the spatial three- vectors and the metric signature is (—+++), so that 

= — Po- K = (A:, k) = (A;, 0,0, A;) is the lightlike momentum of the photon, which 
we choose to be oriented along the z axis. We furthermore assume k > T, which is the 
validity region of the LO and NLO calculations. We will work perturbatively in the strong 
coupling g, meaning that we treat the scale gT (the soft scale) as parametrically smaller 
than the scale T (the hard scale). 

Throughout the paper we will often use light-cone coordinates, which we define as 
p~ = — p^ and p~^ = ^ . This normalization is nonstandard, but we find it convenient 
because dp^dp^ = dp'^dp" , and because we will frequently encounter cases in which p~ = 0, 
in which case p^ = p^ = p~^ with our conventions. The transverse coordinates are written 
as p_L, with modulus p±. 

We finally remark that for convenience we will mostly work in the Keldysh, or r, a , 
basis of the real-time formalism for the computation of Eq. (2.2). The two elements of this 
basis are defined as (pr = (0i+i?!>2)/2, (j)a = 4>i — 4>2, 4> being a generic field and the subscripts 
1 and 2 labeling the time-ordered and anti-time-ordered branches of the Schwinger-Keldysh 
contour respectively. The propagator is a 2 x 2 matrix, where one entry is always zero and 
only one entry depends on the thermal distribution, i.e., 




(2.1) 



a 




(2.2) 




n{p^)){DR-DA) Dr 
Da 




(2.3) 
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where Dn and Da are the retarded and advanced propagators, the plus (minus) sign refers 
to bosons (fermions). n{p^) is the corresponding thermal distribution, either ng{p^) = 
(exp(p^/T) — for bosons or np{p^) = (exp(p''/T) + l)~^ for fermions. We also define the 
spectral function as the difference of the retarded and advanced propagators, p = Dr — Da- 
We will denote the gluon propagator by G and the quark one S. 

We will adopt strict Coulomb gauge throughout. The treatment of soft momenta 
in propagators and vertices requires the use of Hard Thermal Loop (HTL) resummation 
[23]. For convenience we list the Coulomb gauge retarded HTL resummed propagators for 
fermions and gluons in App. A. We will discuss the power-counting rules of the HTL theory 
with fermions in r, a basis in Section 4. 

2.1 Leading-order calculation 

At leading order (and also at NLO) in g the photon production rate arises from diagrams 
where the photon attaches to a single connected quark loop with a number of gluon lines. 
We denote the momenta flowing in the quark lines that attach to one of the photon vertices 
P and K + P - see for example Fig. 1. Then the leading order rate arises from three 
kinematic regimes: 

1. One of the quarks attaching to the photon is on-shell, {K + P)^ ~ g^T"^, and the 
other one is far off-shell, P^ ~ T^. This is the hard 2 -f-)- 2 region. 

2. The photon attaches to one hard on-shell, {P + K)"^ ~ g'^T'^, fermionic line and one 
soft fermionic line with P^ ~ g'^T'^ and P • u ^ gT, where u is the rest-frame of the 
medium. This is the soft 2 -f-)- 2 region, and it is the soft limit of the hard 2 -H- 2 
kinematic region. 

3. The photon line attaches to two fermionic lines which are hard (P • n ~ T and 
{K + P) -ur^ T), nearly coUinear {K ■ P g'^T'^), and nearly on-shell (P^ ^ g^T"^ 
and {P + K)^ ~ g'^T'^); this is the coUinear region. 

The total leading order rate is the sum of these three kinematic regions: 

dr 



dTry 



d^k 



dT^ 



LO 



dr^ 
-I '- 

hard 



soft 



(2.4) 

coll 



The physics of each of these regions will be summarized in the remainder of this section. 

Diagrammatically, the hard 2 -H- 2 region consists of the simple two-loop diagrams 
obtained by adding a single gluon to the one- loop diagram for W^, which is kinematically 
forbidden for a real photon. These diagrams are shown in Fig. 1; their cuts (when all 
lines are hard) correspond to the 2^2 processes qg — t- q'j (Compton) and qq — t- gj 
(annihilation). In Fig. 2 we show an example diagram and its corresponding cut.^ Fig. 3 
shows the processes contributing to the sum of the cut diagrams. 



^In our graphical notation the double line represents particles whose momentum is hard, i.e., 0{T) 
or larger, in at least one component. This thus includes not only hard particles, but also coUinear and 
semi-collinear ones. 
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K + P K + P 




P+Q P+Q 



Figure 1. Two- loop diagrams necessary for tlie evaluation of the 2 -s-> 2 region. The wavy lines 
are photons, curly lines are gluons and plain lines are quarks. No assumption is made yet on the 
scaling (hard, soft, collinear...) of the internal lines. The momentum assignments shown here will 
be used throughout the paper. 




Figure 2. The two- loop diagram on the left corresponds to the square of the amplitude of the 
diagram on the right and on the squares of its crossing. The interference terms arise from the 
two-loop diagram where the gluon is exchanged between the two fermionic lines. Double internal 
lines stand for hard particles, i.e., particles whose momentum is 0(T) in at least one component. 

2 




-I- crossings 



Figure 3. The squared matrix element for 2 o 2 processes. The corresponding annihilation dia- 
grams are obtained by crossing and are not shown. 

The two first (self-energy) diagrams in Fig. 1 receive parametrically equal contributions 
from all logarithmic momentum scales (up to T) of the intermediate virtual quark with 
momentum P. Thus, the naive computation of the diagrams of Fig. 2 leads to a logarithmic 
infrared divergence. However, this IR divergence is regulated by the physics of screening 
and collective plasma excitations, and indeed when the momentum P becomes soft, P ~ 
gT, the self-energy insertion is not anymore a ^-suppressed perturbation in the dispersion 
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relation of the intermediate quark and needs to be resummed. This is the soft region. In 
this kinematic region, the third diagram in Fig. 1 gives a subleading contribution and thus 
the diagrams that contribute to the leading order result are those shown in Fig. 4. These 
diagrams can be most conveniently resummed in Hard Thermal Loop -effective theory, and 
in particular the leading order HTL-resummed diagram is displayed in Fig. 4. The cuts of 
the diagram correspond to conversion processes, where a soft fermion exchange with the 
medium converts a hard quark with momentum K + 0{gT) into a photon with momentum 
K. The computation of the resummed diagram will be discussed in detail in Sec. 4, where 
a sum rule for its evaluation will be introduced. 

The soft and the hard regions are smoothly connected and when gT <^ P <^ T, each 
set of diagrams gives an equivalent and correct description. In a practical calculation one 
introduces a momentum cutoff for the integrals gT <^ fJf^ <^ T dividing the two regions. 
Both contributions then individually depend on fi^, but this dependence cancels exactly 
in the sum of the two terms yielding a logarithm of the temperature over the asymptotic 
mass moo ~ gT of the quark. 



Wk 



soft 



+ 




+ 



Figure 4. The diagram contributing to the soft region at leading order. The double lines are the 
hard lines whereas the dotted single line represent bare soft propagators. The single plain line is 
the HTL-resummed soft propagator. The momentum assignments are given in Fig. 1. 

In [13, 24, 25] it was pointed out that there is another phase space region that con- 
tributes at LO besides the hard and soft regions: the coUinear region. In this region 

• the momenta of the quark lines are hard, nearly on shell and collinear to each other, 
i.e., T, p± gT, and p~ ~ g'^T such that ~ g'^T, and 

• the momentum Q of the gluon is spacelike and soft with g+ ~ g_L ~ gT and q~ ~ g'^T, 
so that the kinematics of the quarks are unaffected by the gluon. 

These constraints force the angles between the quarks and the photon to be small ~ g, 
and therefore the cuts of the third diagram in Fig. 1 correspond to the hremsstrahlung and 
pair annihilation processes shown in Fig. 5. In this case the intermediate virtual quark is 
almost on shell and thus has a long lifetime of order ~ ^/g^T which is parametrically of the 
same order as the small angle scattering rate in the plasma l/F ~ 1/g'^T. Hence during 
the formation of the photon, there is an 0(1) probability that the intermediate quark 
will undergo one or several additional soft scatterings with the constituents, and processes 
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Figure 5. CoUincar diagrams. In the first case, called the bremsstrahlung diagram, the angle 
between the emitted photon and the outgoing emitting fermion is of order g. In the second case, 
called the pair annihilation diagram, it is the angle between the annihilating quark and antiquark 
that is of order g. The diagrams where the gluon is attached to the other fermionic line are not 
shown. In both graphs the gluon is soft and is scattering off the hard quarks and gluons of the 
plasma as indicated by the crosses, i.e. it is an HTL gluon in the Landau cut. 



involving multiple soft scatterings are not suppressed by powers of g. These multiple 
scatterings lead to destructive interference, that is known as the Landau- Pomeranchuk- 
Migdal (LPM) effect that leads to an 0{1) suppression of the collinear rate. 

In terms of the two-point function these processes correspond to diagrams with the two 
nearly collinear fermion lines connected with arbitrary number of soft spacelike gluons with 
same kinematics as Q. In [14, 15] Arnold, Moore and Yaffe (AMY) showed that it is only 
the ladder-type diagrams shown in Fig. 6 that contribute to a leading order calculation; the 
factors of g arising from additional vertices are canceled by near on-shell propagators and 
large statistical factors arising from the gluonic propagators. The near on-shellness of the 
quark lines makes the diagrams sensitive to the thermal mass ~ g^T"^ and the thermal 
width r ~ g^T of the quark lines, which need to be consistently resummed. Furthermore 
AMY showed how these diagrams can be resummed in terms of a Schrodinger equation 
type differential equation, and they obtained the complete leading-order results in [15]. In 
Sec. 3 we will discuss in detail this equation in the context of the treatment of its NLO 
corrections. 



Wk 



coll 




=Re 
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Figure 6. The uncrossed ladder diagrams that need to be resummed to account for the LPM effect 
in the collinear region. The cut shown here corresponds to the interference term on the right-hand 
side. The rungs on the l.h.s. are HTL gluons in the Landau cut. On the r.h.s., the crosses at the 
lower end of the gluons represent the hard scattering centers, either gluons or fermions. 
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The leading order result arising from the hard, soft and collinear regions can be sum- 
marized as [15] 



dV. 



d^k 
dV^ 



d^k 



hard 



soft 



coll 



A{k) 

(27r)3 

A{k) 

(27r)3 

A{k) 

(27r)3 



In 



In 



+ Chard 



T 



'-'coll 



or 



dV^ 



d^k 



LO 



A{k) 
(27r)3 



In 



rrir 



+ Chard ( ^ j + C^ou [ ^> 



(2.5) 
(2.6) 
(2.7) 

(2.8) 



where A{k) is the leading-log coefficient of the photon production rate 



Aik) = 2aEM 



m 



AaEMnF{k)g^T^ 



for QCD with uds quarks 



k 'Z_^— ^ 

(2.9) 

Here d^ is the dimension of the quark's representation (d^ = Nc in the fundamental 
representation of SU{Nc)), Qs is its abelian charge, and the sum runs over the number of 
light fermions flavors, Nj. The parameter k is the square of the ratio of = g^T'^Cf./A to 
the Debye mass m\ = g^T"^ {C a + N jTr) /2> at leading order in g. k encodes the dependence 
on the number of colors and light flavors in the plasma, i.e. 



3C« 



rrij 



4{Ca + TRNf) 



6 + N 



f 



for QCD 



(2.10) 



where Cr and Ca are the quadratic Casimirs of the representations of the quarks and gluons 
respectively, and Tr is the index of the representation of the quarks.^ The definition of 



is discussed in Appendix B.4. 



The functions Chard(^/r) and C^^(A;/T, k) describe the momentum dependence of the 
hard and collinear regions and have to be obtained numerically, the former by integrating 
the matrix elements for the hard processes folded over the thermal distributions and the 
latter by solving the integral equation for collinear processes. For further convenience we 
list the parametrization of these functions as given in [15] for a Nc = 3 QCD plasma with 
Nf flavors: 

Chard (^) - I In + 0.04l| - 0.3615 + l.Ole-'-'^'^/^ , 0.2 < ^ (2.11) 



c 



TO?1 [ f^'^ 



0.316 ln(12. 18 + r/A;) 



(k/T) 



3/2 



+ 



0.0768fc/T 
^1 + A:/(16.27'fy 



0.2 < - < 50. 



(2.12) 



^ Cr = (N^ — l)/2Nc and Tr ~ 1/2 for quarks in the fundamental representation, and Ta ~ Ca ~ Nc 
for the adjoint representation of SU{Nc). 
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We conclude this overview of the leading-order calculation by noting that the momen- 
tum integration regions that contribute here are best identified by their scaling in terms 
of P. In Fig. 7 we map these momentum regions in the {p'^tPa_) plane. The scaling of p~ 
can be obtained from momentum conservation. 




Figure 7. Momentum integration regions in the {p'^ ,p\^ plane contributing to the leading-order 
calculation. The ^i^^ label indicates a LO cancellation of UV/IR log divergences between the soft 
and hard regions respectively. 



2.2 Next-to-leading order corrections 



At next to leading order, the full result is a sum of the leading order rate and its 0{g) 
correction 



LO+NLO 



+ 



LO 



d5T^ 



(2.13) 



As in the leading order calculation, the NLO rate arises from distinct kinematic regions 
and the NLO correction can be parametrized as 



d5T^ 



d6T^ 



soft 



d^k 



+ 



d5r^ 



coll 



d^k 



(2.14) 



semi— coll 



The soft and collinear regions are the same kinematic regions as in the leading-order calcu- 
lation, while the semi- collinear region is an additional kinematic region whose contribution 
starts at NLO. 
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In the hard region, corrections come about by adding an extra loop to the diagrams 
shown in Fig. 1. As long as momenta stay hard, we need not worry about these corrections, 
which are suppressed by a factor of relative to the leading-order result. 

However, both the soft propagator and the location of the quasiparticle pole at hard 
momenta have 0{g) corrections, and thus diagrams that are sensitive to these quantities 
may also receive 0{g) corrections. Indeed, in the collinear case the individual rungs in 
the ladder resummation are soft gluons, and the collinear quarks are near the quasiparticle 
pole. Consequently, the parameters of the collinear integral equation that record the soft 
gluon scattering rate and the quasiparticle masses, C(g^) and m^, are modified at NLO. 
The structure of the NLO correction arising from the collinear region is then 

d5T. 



d5T^ 



d6T^ 



coll d^k 



+ 



(2.15) 

5C 



where the first term is due to a 0{g) shift in m^, and the second term arises from a 
one-loop 0{g) correction to the soft scattering kernel, C{q±). In Section 3 we will discuss 
these corrections in detail, as well as two O (g) subregions where the collinear phase space 
starts to overlap with the soft and semi-collinear regions, requiring subtractions. 

In the soft region, the addition of an extra soft gluon to the diagram in Fig. 4 results 
in the diagrams shown in Fig. 8, which represent an O (g) correction to the conversion 
process. In particular, wherever a gluon ends on a soft fermion line, all momenta fiowing 




Figure 8. Diagrams contributing to the NLO fully soft rate. The black blobs are bare+HTL 
vertices, plain lines and gluons are soft. We call these four diagrams, from left to right, the soft-soft 
self-energy, the tadpole, the hard-soft self-energy and the cat eye. The momentum assignments are 
given in Fig. 1. 

in that quark-gluon vertex are of order gT. This causes the bare and HTL vertices to be 
of the same order, requiring the inclusion of the HTL vertex, as shown in the first and 
last diagrams in Fig. 8. Furthermore, the two-quark, two-gluon HTL gives rise to a new 
topology, the second diagram in that figure. 

The complicated analytic structure of the HTL vertices and propagators, with their 
branch cuts and imaginary parts, as well as the non-trivial functional dependence on the 
momenta, would in principle make the calculation of the diagrams in Fig. 8 technically 
intricate and only amenable to a multi-dimensional numerical integration. However, in 
Sec. 4 we develop a set of sum rules using the analytic properties of these amplitudes, 
which are in turn related to causality. These sum rules, as we shall show, simplify the 
calculation dramatically leading to an analytical result. 

When evaluating these soft diagrams we must correct the LO treatments of the soft 
region to avoid double counting. Note that the first diagram in Fig. 8 is the soft limit 
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of the HTL self-energy already included in the soft-LO calculation, see Fig. 4. The HTL 
self-energy used at LO includes an integral over the hard thermal loop momentum Q 
which extends down to zero, with 0{g) of the contribution arising from the soft region of 
integration where the HTL approximation is no longer valid. Thus, while the first diagram 
in Fig. 8 has already been included at LO, the LO treatment is incorrect at NLO. To avoid 
double counting and to correct this mistreatment we have to subtract this soft-loop part of 
the HTL calculation. We do this by subtracting a counterterm diagram shown in Figure 
9. 




Figure 9. Mistreated 0{g) part of the leading order soft calculation. The dotted lines indicate 
bare soft propagators. 

Similarly, in the calculation of the leading order collinear rate, an 0{g) part of the p"*" 
integration arises from the kinematic region where p"*" is soft and overlaps with the soft 
region. For example, the LO treatment of the collinear bremsstrahlung diagram shown 
Fig. 5 integrates over the momentum fraction of the final state quark, and is incorrect when 
this momentum fraction is 0{g). This region is correctly dealt with by the soft graphs of 
Fig. 8, and the mistreated LO collinear contribution must be subsequently subtracted. The 
structure of the soft region at NLO is then 







d5T^ 


diags. 


dV^ 


subtr. 




soft 


d^k 


soft 


d^k 


soft 



(2.16) 



where the first term arises from the difference of the diagrams in Fig. 8 and Fig. 9 and 
the second is the soft part of the leading order collinear rate. The complete treatment of 
the soft region, including the diagrams of Fig. 8 and all necessary subtractions, is given in 
Section 4. 

Finally, a third region contributes to the NLO rate, corresponding to the uncharted, 
wedge-shaped area in Fig. 7 between the hard and collinear regions where Q ~ gT is soft 
and the scaling of P obeys ^ T, p~ ^ gT and p± ~ ^fgT ^ so that ~ gT"^ ^ P ■ K 
— this is the semi- collinear region. This region is closer to the mass shell and closer to 
collinearity with K than the hard region (where P"^ P ■ K), but farther from the 

mass shell and less collinear than the collinear region (where P^ ~ g'^T^ ~ P ■ K).^ The 
physical processes in this region are characterized by the sign of Q^. For timelike gluon 

^Naive power-counting arguments suggest that this region, witli the exchange of a soft Q ~ gT gluon, 
should actually be leading order. But there is a cancellation between the diagrams shown in Fig. 1. This 
cancellation, pointed out in footnote 5 of Ref. [14], renders this phase space region 0{g) with respect to the 
leading contributions. 
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momenta the semi-collinear contribution can be seen as a correction to the hard region 
from the emission/absorption of soft gluons on their massive, quasiparticle plasmon poles. 
In the spaceUke region, on the other hand, one can interpret this contribution as a less 
collinear emission, the angle being now yjg instead of g. This in turn implies a formation 
time of the order of \/{gT), much shorter than the inverse scattering rate. The LPM effect 
is then not relevant in this region. 

The leading order treatment of the collinear region also mistreats the semi-collinear 
region by an 0{g) amount, and the part of the leading order collinear rate arising from the 
semi-collinear region again needs to be subtracted. Also, the leading order hard calcula- 
tion receives an 0{g) mistreated contribution from the semi-collinear regions so that the 
structure of the semi-collinear correction becomes 



d5T^ 



semi— coll 



diaes. 



semi— coll 



dT^ 



d^k 



coll. subtr. 



semi— coll 



dr.. 



d^k 



hard subtr. 



semi- coll 



(2.17) 



The complete treatment of the semi-collinear region is given in Section 5. 

We conclude this section by mapping the NLO integration regions in the plane 
in Fig. 10, in analogy the LO map in Fig. 7. 




Figure 10. Momentum integration regions in the {p^,p±) plane contributing to the LO and NLO 
calculations. The yellow and pink bands indicate the LO and NLO integration regions respectively. 
The blue bands indicate an overlap region where the LO contribution must be subtracted to avoid 
double counting. The hard region does not contribute at NLO. Finally, the label anticipates 

a cancellation of UV/IR log divergences between the soft and semi-collinear regions, while the 
label indicates a similar LO cancellation between the soft and hard regions. 
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3 The collinear region 



The evaluation of the colhnear region at leading order requires the resummation of an 
infinite number of soft gluon exchanges through an integral equation. Such an equation 
was derived by Arnold, Moore, and Yaffe [14, 15] and gives rise to a LO contribution to 
the photon production rate of 



d^k 



coll 



Ajk) 
(27r)3 



dp^ 



(p+)2(p++A:)2 



21 



np{k) 



g^C^T^J (27r)2 

d'q± 



2px = i6E f(px) + 

kjpl + m^) 
2p+{k+p+) 



(2vr) 



Re2p_L • f(p_L,p"^, A;) , 
2 C(gx)[f(p±)-f(p±+q±) 



6E 



(3.1) 
(3.2) 
(3.3) 



Here 5E = -Bp+k ~ Ep — k is the eikonalized energy difference between having a quark of 
energy corresponding to a momentum oi k+p and having a quark of energy corresponding 
to momentum p with a photon of energy k. C{q±) is the differential soft scattering rate, 
which at leading order reads [14, 26] 



dq^dqz 
(2vr)2 



27r5(g°-g,)G-(Q)<^^ = 5'C7,r 



rrir 



qliql+ml) 



(3.4) 



where 



1,0,0, 1) and GJ^^ is the cut HTL gluon propagator. 
The physical interpretation of this expression is as follows. Photon production involves 
a current operator insertion in the amplitude, followed by time evolution, and then a current 
insertion in the conjugate amplitude. At times between the current insertions, the density 
matrix contains an off-diagonal term with a quark with momentum p~^ , p± and a photon 
of momentum k in the amplitude, but a quark of momentum (p'^+k), p± in the conjugate 
amplitude. The size for this entry in the density matrix is f(p_L) and Eq. (3.2) is the time- 
integrated evolution equation for this amplitude; 5E represents the phase accumulation 
because the state in the amplitude and conjugate amplitude do not have the same energy, 
while the C{q±) term describes the effect of scattering processes on the evolution of the 
density matrix. The second line in Eq. (3.1) describes the overlap of the current operator 
on this density matrix element. The term in square brackets in the first line of Eq. (3.1) is 
the DGLAP splitting kernel. 

The NLO corrections to this leading order calculation arise from an 0{g) correction 
to thermal mass and to the differential soft scattering rate. These contributions, like the 
LO contribution, need to be solved for numerically. The most convenient way to do so 
is to Fourier-transform Eq. (3.2) to impact parameter space, where the integral equation 
turns into a differential equation with mixed boundary conditions. This will be discussed 
in subsection 3.2. First we find the behavior of Eq. (3.1), Eq. (3.2) at small p"*" and at 
large pj_, which we need as counterterms in Eqs. (2.16) and (2.17). We can find these 
analytically by perturbing Eq. (3.2) in 6E. 
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3.1 Leading order subtractions: |scmi-"oif' -'^7 1 soft*'^' 

In the leading order calculation of the collinear rate in Eq. (3.1), the integral over 
extends artificially all the way down to < gT; and the p± integral extends up to 
p± ^ gT, contrary to the definition of the collinear kinematic region. In these regions, 
the integrals start to probe the soft and semi-collinear regions, respectively. There, the 
assumptions made in arriving at Eq. (3.1) and Eq. (3.2) start to break down, so that a 
more detailed analysis is needed. To find the leading order contribution we simply extend 
the p"*" integral to zero and the p± integral to infinity; then when we handle the p'^ ~ gT 
and ~ gT^ regions more carefully, we will subtract the (incorrect) amount already 
included in this way in the leading-order calculation, via the subtractions of Eqs. (2.16) 
and (2.17). So let us evaluate Eq. (3.1), Eq. (3.2) in these regions. 

In the soft and semi-collinear regions, one finds that 6E ~ gT, and therefore in these 
regions the first 6E term in Eq. (3.2) is larger than the second C{q±) term. Hence, the 
evaluation of Eq. (3.1) and Eq. (3.2) can be simplified by working to first order in C{qj_). 
Physically this is because, if the opening angle becomes large 0{yfg) or if an external quark 
becomes soft and therefore scatters at a large angle, then the formation time becomes ~ 
1/gT, parametrically shorter than the mean time between scatterings ~ 1/g'^T. Therefore 
the LPM effect becomes subleading and the rate is determined by the single scattering 
rate. 

In order to solve Eq. (3.1) by perturbing in write f = fi-|-f2-|-. . . with fi oc ((^i?)~^, 

f2 oc {8E)~'^ etc. and evaluate iteratively: 

2p^ = i5i?fi(p^) ^ fi = -|B^, (3.5) 

which gives zero in Eq. (3.1) because it is imaginary. Substituting fi into the collision term 
to determine f2, we get 

= i<5i?(px)f2(p±) + j ^C(gx) [fi(p±) - fi(p± + q±)] , 

^^^P-) = jEi^) J (^^(^-) [sm^) - 5^(p.+q.)J • ^'-'^ 

Since this term is real, it contributes to Eq. (3.1). The next term £3 is again imaginary, 
and £4 is 0{g) in the regions of interest and therefore negligible. Integrating over p± as in 
Eq. (3.1) and symmetrizing the resulting expression with respect to p±, {p±+q±), we find^ 

f d'^P^ n CI \ ^fd^3±d^r/ ^/^ P± + q± V ^Q7\ 



''This symmetrization is natural, indeed necessary, from the point of view of the momentum labeling 
which we introduced on the diagrams shown in Fig. 1. Namely, the l/{SE{p±))^ term corresponds to 
the first diagram there (self-energy on the lower line), the l/{5E{p±+qj_))'^ term is the next diagram 
(self-energy on the top line), and the cross-terms are the cat eye diagram. In terms of Fig. 5, for the 
bremsstrahlung diagram the l/{5E{p±))'^ term is the diagram shown, the l/(5iJ(p_L-|-qx))^ term is the 
square of the diagram where the gluon attaches after the photon, and the cross-term is the interference 
term between these diagrams. Not performing the symmetrization corresponds to evaluating one of the 
diagrams with different momentum assignments on the external lines. 
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We now turn to the application of this equation to the specifics of the soft and semi-colhnear 
hmits. 



3.1.1 The soft fermion, collinear contribution 

We now consider Eq. (3.1) in the region where either or p"*" + is small. We will 
show that the integrand goes over to a constant, so 0{g) of the contribution to arises 
from the region where is 0{g). This ^'''-independent behavior will turn into a linearly 
divergent subtractive counterterm when we evaluate the soft region. 

First consider Eq. (3.1) in the regime where is formally 0{T) but soft, p~^ <^T and 
p'^ ^ k. In this case we can approximate Eq. (3.1) as 



(Pk 



- ^^!:]. .^rfP^TzL^ I |^Re2p^.f(p^,p+). (3.8) 

soft 



(2^)3 J ^ {p+y 2 g^C^T^ J (27r 



Upon plugging the soft-p''' limit of 5E, i.e. 

2p 



6E = , (3.9) 



into Eq. (3.7) we obtain 

J (2vr)2 J (27r)4 \^p^+r7i^ (p±+q±)^ + 

(3.10) 

In terms of p"*" scaling, we see that Eq. (3.2) gives {p'^)'^ times a ^^-independent 
function. This cancels the (p^)~^ in the integrand in Eq. (3.8), so indeed the integrand in 
Eq. (3.8) is independent of p~^ at small p"^. Since p~^ ~ gT represents 0{g) of the phase 
space of p"*" values available, this region therefore represents an 0{g) fraction of the photon 
production rate, as claimed. 

The region where p'^+k is soft gives an identical contribution. Inserting 2f2 into 
Eq. (3.8) we then get 



d^k 



A{k) /■+^+ ^ + 8 /■ dPp^cPq^ ml 



dp — 



(2vr)3y_^+ - Tj (2vr)4 g^gi+m^) 

p± p± + qj 



.PI + (p±+q±)^ + 

where we introduced a regulator gT <^ /i'^ ^ T for the linear divergence. 



(3.11) 



3.1.2 The semi-coUinear fermion, collinear contribution 

The semi-collinear region represents another 0{g) contribution to the integral in Eq. (3.1). 
As in the previous case, the approximations that lead to that equation are no longer valid 
when P becomes semi-collinear (p± — t- y/gT, p^ ~ gT). This limit is then incorrectly 
described by Eq. (3.1) and, as in the previous subsection, we need to derive its limit in 
order to subtract it from the semi-collinear region, where this momentum scaling will be 
correctly treated. 
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We can again use Eq. (3.7), but now there is an additional simplification; ^ 
and I ^ I I . Therefore we can drop and work to lowest order in , which is 



^2p± ■ f2(p±) = 2 / ^7T3^C(g± 



(2vr)^ 



^1 



(27r)4 ^^-'(<5ii;(p^))2 

4(p+)2(p+ + A:)2 f <fq^ ^ 



(2vr)2 



(2vr) 



(?iC(gx). (3.12) 



When plugged in Eq. (3.1), this yields 



d5T^ 



coll subtr. 



semi— coll 



2^ / dv^ 
^(27r)3 ' "'^ 

1 



(p+)2 + (p+ + ky 



np{k) 



g^C^T^J (27r)2 



d^L A{p+f{p+ + kf f d^q± 2 



(2vr)' 



■(?lC(gx). (3.13) 



The p± integration is power pj_ divergent and the qj_ integral is log UV divergent. This is 
not surprising, since this expression was obtained based on qj_ ^ 



3.2 NLO corrections to the coUinear regime: 6r^\sm and ^F^l^c" 

Even at leading order Eq. (3.2) has to be solved numerically in order to get the collinear 
contribution. The most convenient way to do so is by Fourier transforming p± and qj_ 
into impact-parameter variables, as first proposed in [27]. The advantages are, first, that 
the convolution over the collision kernel C{q±) becomes a product, turning an integral 
equation into a differential equation; second, that the source on the left-hand side becomes 
a boundary condition at b = 0; and third, that the desired final integral, Eq. (3.1), becomes 
a boundary value of the ODE solution. Specifically, defining 



we have 



and Eq. (3.2) becomes 



Re 



f(b) 



(2vr)2 



/ 



(2vr)2 



f(q±) 



2p± • f(p±) = Im(2Vb • f(6)) , 



ik 



2p+{k + p-^ 



wi 



th^ 



(2vr)2 



1 



vnf(b)+c'(6)f(b) 



Ciq±) 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



^Note that C'(6) is not the Fourier transform of C{q±), but rather the difference between the Fourier 
transform at zero b and at finite b, which is better behaved (in particular, not sensitive to the divergent total 
cross-section). Alternately, we can redefine C{q±) to have a negative delta function at q± = 0, normalized 
so that its integral J d^q±C{q±_) vanishes, in which case Eq. (3.2) does not need the first term in square 
brackets, and C'{b) is minus the Fourier transform of C{q±). 
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In the collinear regime the 0{g) corrections enter then in two places: both the effective 
thermal mass squared and the collision kernel C{q±) get 0{g) corrections which modify 
Eq. (3.16), 



m 



oo,LO+NLO 



+ 6m^ 



LO+NLO 



(6) = C'{h) + 5C'{h). 



(3.18) 
(3.19) 



The computation of the NLO thermal mass from [22] is rederived in Appendix B.4. The 
NLO collision kernel is computed in [21] in momentum space; in Appendix C we perform 
the Fourier transformation into impact parameter space. 

Eq. (3.16) is then solved perturbatively, by treating f(b) formally as an expansion in 
powers of 5m^, JC; f (b) = fo(b) + fi(b) + . . ., and expanding to first order. The zero-order 
expression is just Eq. (3.16), while at the linear order the expression reads 

ik 







ik 



vn+c'(6) fi(b) + 



+ 5C'(6) fo(b), (3.20) 



2p+{k+p+)\ ~ "J wy XV y > \2p+{k+p+) 

where the leading order solution fo(b) acts as a source term in the diff'erential equation for 
fi(b). 

When evaluating Eq. (3.16) one must deal with mixed boundary conditions; the func- 
tion f(b) must decay to zero as b — )• oo (one boundary condition), and it must yield the 
correct normalization 

72._,u^_4p+(fc+P+)^,2, 



V^fo(b) 



k 



(3.21) 



at zero. This is done by evolving the differential equation starting at large b, with starting 
data which ensure that the solution will decay away as b — t- oo but with an arbitrary 
normalization. One then solves the differential equation going in towards the origin, gen- 
erally resulting in a mis-scaled solution. But this solution can be multiplied by a complex 
constant so that Eq. (3.21) is satisfied. Similarly, when solving Eq. (3.20) for fi(b), the 
boundary condition that fi(b) should decay at large b is not enough to fix the solution 
completely; so one generically gets a solution which is a mixture of the solution to Eq. (3.20) 
with correct boundary condition linib-j-o V^fi(b) = 0, plus a multiple of the homogeneous 
solution (that is, the solution of Eq. (3.20) at fo = 0) with the wrong boundary condition 
at zero. But the homogeneous solution is proportional to fo, which is known; so it can be 
subtracted to obtain the solution with correct boundary condition. 

We solved Eq. (3.20) as a function of k/T. The results can be parametrized as follows: 



dST^ 



d?k 
d5T^ 



d^k 

d5T^ 
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(5m 
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d^k 

A{k) 

(27r)3 
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(27r)3 
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d5T^ 
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^^oo r<5m 
7\ — 



5C 
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nir 
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2m£, 



(3.22) 
(3.23) 
(3.24) 

(3.25) 
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Figure 11. Two NLO functions, ((5m^/m^)Cf^i(fc/r) and {g^CAT/mr,)Ci^y^{k/T), which pa- 
rameterize the changes in the the colhnear emission rate due to the NLO quasi-particle masses and 
colhsion kernel respectively - see Eq. (3.22). The full LO+NLO coUinear emission function is a sum 
these two corrections and the leading order result, C]^2\ + '^C'coii- The curves are for Nc = Nf = 3 
and Us — 0.30 

The fitting functions read 




- 0.3664/x - 0.08478 - 0.0799 log(x) 



+0.0315 X - 0.0050 X log(x) - 0.0681 {\og{x)f^ , (3.26) 

Cf^j(rE) = (-0.7207 - 0.82365k + 3.9865k2)/x 
+(0.7056 + 0.0998(5«; - 1.1865^^) 
+(-0.8309 - 0.2610(5fi; + 2.247.5^2) log(a;) 
+(0.12305 - 0.01085k - 0.2871.5^2) x 
+(-0.01777 + 0.00148<5k + QMU5n^)x\og{x) 
+(0.2804 - 0.0369(5k - 0.2375(5k2) \og{xf 

+(-0.0702 + 0.004405k + 0.11495k2) \og{xf , (3.27) 

with 5k = K — 2/9. The fitting functions have a relative error smaller than 2% for Nf = 3 
QCD (k = 2/9) in the momentum range 0.5 < a; < 70. In the range 0.15 < k < 0.35 the 
relative error is less than 5%. 

We will present most of the numerical results, for different values of the parameters 
such as the coupling, in the exposition of the final results in Sec. 6. Here we just show in 
Fig. 11 the size of the mass correction, [5m^/m^]C^™j(/c/T, k), and the collision kernel 
correction, [(72CAT/m£,]C^^[[(A;/T, k), relative to the the LO collinear result, (^^^^[(/c/T). 
The NLO correction 6Cco\\{k/T) is an 0(100%) correction for most of the considered range. 
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Wr 



Wa 




Figure 12. Leading-order diagrams in the r, a basis. The lines witli arrows indicate retarded ra 
propagators with the arrow flow from a to r, while lines without arrows indicate rr propagators - 
see the text for further discussion. The first two diagrams contribute to Wr, with either the bottom 
or the top propagator retarded; the last two contribute to Wa- The diagrams with cut soft lines 
(second and fourth diagrams) arc suppressed by the small statistical function on the cut soft line. 

4 The soft region 

To introduce the NLO calculation, we begin by reproducing the soft-momentum part of the 
leading-order calculation introducing the notation. We perform the leading order calcula- 
tion using novel sum rule technology, which admits a generalization to the NLO calculation. 

4.1 Leading-order evaluation and introduction to the fermionic sum rules 

The most straightforward approach to the leading-order calculation is to evaluate the 
trace of the Wightman correlator W^{K) = g^^'^W^^{K) in Eq. (2.1) directly in the 1,2 
basis. However, it turns out that the NLO calculation is much simpler to perform in the 
r, a basis, so we will instead use the r, a basis also at leading order. The object that is 
most conveniently calculated in the r, a basis is the retarded (advanced) correlator, which 
is related to the backward Wightman function via the KMS relation 



At leading order, Wr{K) and Wa{K) each arise from two r, a assignments of the 
one loop diagram, shown in Figure 12. Our graphical notation for the r, o assignments 
follows the one in [28]: we draw outgoing arrows for a fields at vertices and incoming 
ones for r fields. For the rr propagator, rather than drawing two arrows pointing in 
opposite directions, we omit to draw them; this should cause no confusion. The double 
lines refer to hard {K + P) T propagators whereas the single lines refer to soft P ~ gT 
HTL propagators. The different r, a assignments of resummed HTL propagators are easily 
obtained from the retarded ones listed in App. A by using Eq. (2.3). 

The cut (rr) soft line carries a factor of — n^(p'^)-|-l/2 <^ 1 and is therefore suppressed;^ 
we may therefore drop the two diagrams containing cut soft lines. Summing the other 
diagrams gives the difference between retarded and advanced propagators on the soft line. 
This difference is the spectral function p = Sr — Sa- Also the cut line can be expressed in 
terms of the spectral function; Srr{K + P) = (| — np{kP +p^))p{K + P). We approximate 

®It is also odd in and will be even more suppressed when averaging over p+ — >■ — p"*". 



{K) = 2ns (A;) l^iWR{K) 
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this statistical function as ^ — np[k^ + p^) — \ — np{k) and use the identity 

neik) {I - 2np{k)) = np{k) . (4.2) 

Bringing everything together, 

^<(^) ^ ^ gyd^Mk) I |p ^ j^^^^^ ^ ^ ^^^^ 

If we had evaluated directly without going to the r, a basis we would have written this 
down immediately. 

Now we evaluate Eq. (4.3) expanding in P ~ gT K ^ T. This expansion will 
enforce eikonality on the hard line, which is essential for the sum rules described below. It 
is convenient to write each propagator in terms of its components of positive and negative 
chirality-to-helicity ratio: 

p{P) = h%p^{P) + h-p-{P), (^-^^ 

with p = p/p. For the hard line, we use an eikonal approximation 

p-^{P+K)c^2T:6{vk-P) = 2T:6{p-), p-(P+/^) = 2^(5(p° + A; + |p+k|) ~ , (4.5) 

and thus the hard line is a function of p~ only. Using this delta function simplifies the 
traces, which we expand in small P: 

T. = 2 (l . ) ^ 2| ± p^M + o (g) , (4.0) 

We insert the leading-order piece of this trace into Eq. (4.3), finding 



W<{K) = 2Y,qydnnp{k) j 



(27r)2 



(4.7) 

The factor of two in front takes care of the kinematic region where K -\- P \s soft and P is 
hard, which gives the same result after a shift of integration variables. 

Next consider the p"*" integration in Eq. (4.7). We will perform the integral using 
analyticity methods similar to those discussed in App. B. The key is that Eq. (4.7) in- 
volves p{p~ iP'^ ,Pa.) = Sr(j>~ ,p'^ ,p±) — Sa{p~ ,P'^ ,p±)- Due to causality, the retarded and 
advanced functions are analytic in any timelike or null momentum variable in the upper 
and lower half planes respectively, generalizing the familiar analyticity properties of these 
functions in p^ [21]. In particular, Sji{p~ ,p'^ ,p±) is analytic in the upper half of the com- 
plex p~^ plane, while holding p~ and p± fixed. We are therefore free to deform the p~^ 
integration contour: instead of integrating just above and below the real axis for Sji and 
5^1 respectively, we integrate along an arc at large p"*" where gT <^p^ <^T ~ see Figure 13. 
Along these arcs the integrand has a remarkably simple behavior, which can be obtained 
by expanding the HTL propagator listed in App. A, yielding 

^2 



1 - S+{p+,p)+[l + ^ ) S-{p+,p) 



O ( 7^] ■ (4.8) 
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Figure 13. Integration contour in the complex p+ integration, and the deformation we use to 
render p+ ^ gT. Wr runs below the real axis and Wa above. This happens because the letters 
R and A refer to the causal prescription with respect to K and our momentum assignments imply 
that when W is retarded in K it is advanced in p'^ . 

Integrating along the arcs at positive and negative Im(p"'") for the retarded and advanced 
contributions one then obtains 

W<{K) = 2j2qyd,nAk) [ ^^^V" (4-9) 

A numerical integration of Eq. (4.7) agrees perfectly with this expression. This result was 
also recently obtained in [29]. 

To summarize, the reason why it is possible to deform the integration can be 
understood diagrammatically from Fig. 12. When evaluating Wr and Wa (the first and 
third diagrams in Fig. 12), the soft fermionic lines are either fully retarded or fully advanced, 
since soft rr fermionic propagators are suppressed. This is seen from the flow of arrows on 
soft fermionic lines, which indicates the r, a assignments. Thus, since the cut hard line is 
eikonal and is only a function of p~ and not p"*" , the p'^ integration is over either a fully 
retarded or fully advanced function and can be deformed away from the real axis. 

It is worth noting that, while each of the two components {l=fp'^/p)S^{p~^,p) presents 
poles separately at p^ = ±ip± and branch cuts in {—ioo,—ip±), {ip±,ioo), their sum is 
analytic for Im{p~^) ^ 0. This must be so as our sum rule is based on causality; but the 
chirality-to-helicity decomposition is not Lorentz covariant, so the individual terms need 
not respect causality. This can be seen from the properties 

S+{-P) = S~*{P) , 5-(-P) = 5+*(P), (4.10) 

from which one can also see that the sum of the two components is covariant and does 
respect causality. 
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To constrain the integral to the soft region only, the p+ integral should be cut off at a 
finite momentum scale ^ gT. However, if we cut off the integral at a finite fi~^ ^ gT, 
the 0{(j)'^)~'^) subleading corrections to Eq. (4.8) give rise to {gT/ fi'^)"'~^ suppressed 
corrections. But these will be canceled by corrections which will arise when we perform 
the calculation of the region above /x"'", since the total result should be independent. 
Therefore we need not compute them. 

The (Pp± integration in Eq. (4.9) should be cut off at some large momentum gT <C 
<C T, where it should match with the contribution from the hard region. The explicit 
expression reads 



which agrees with the original calculations of the LO soft region in [11, 12]7 

The manipulations made in arriving from Eq. (4.3) to Eq. (4.7) are valid up to NNLO 
corrections, and we do not need revisit them in the NLO computation. While the approxi- 
mations we have made to statistical functions have 0{g) corrections, they give rise to odd 
integrands in and would hence give a vanishing integral when plugged in Eq. (4.7) since 
the term in square brackets is even as given by Eq. (4.10). In a similar way one can show 
that the order-^ correction from the trace, i.e., the order p\/{pk) term in Eq. (4.6), results 
in an odd integration. Finally, we can consider the 0{g) correction to the dispersion rela- 
tion of the hard line, which changes 6{p~) to 8{p~ — {p\ + m^)/(2fe)). But the difference 
between these delta functions again yields an odd integrand at the 0{g) level. Eq. (4.7) is 
then free of 0{g) corrections. 

However, the soft HTL fermion propagator is resummed in the hard self-energy inser- 
tions as shown in Figure 4. Whenever the momentum flowing inside these internal loops 
becomes 0{gT) soft, the approximations made in the computation of the HTL propagator 
fail. The diagrams where exactly one of the internal loops becomes soft represent a mis- 
treated relative 0{g) contribution which needs to be subtracted in the NLO calculation. 
We will return to this contribution in Sec. 4.3.2. 

4.2 The structure of the soft NLO corrections; a quick derivation 

We saw in Sec. 2.2 that the leading-order diagram of Fig. 4 receives ordex-g corrections from 
four diagrams, shown in Fig. 8. These correspond to four different ways to add a soft gluon 
to a LO diagram. In these diagrams HTL corrections appear on all gluon propagators, 
all soft fermionic propagators, and on vertices in three of the four diagrams. We will 
compute each diagram in detail in the following subsections. But here we will present a 
not-quite-rigorous argument which establishes what the sum of the diagrams must yield. 

^The calculation of [12] used a different regularization, cutting off the d?'p/{2-K)'^ integral at p — /i. For 
any UV log-divergent function in three dimensions that for p > /i is approximated by its asymptotic behavior 
l/p^, the difference between our cylindrical regularization {p± < fi, -co < p+ < cxd) and their spherical one 
is (1 — ln(2))/(27r^). Indeed, by inspecting Eq. (17) in [12] and fixing the overall normalization one sees that 
their result is A{k){ln{jj,/moo) + ln(2) — 1) (their numerical term on the second line is —0.31 « ln(2) — 1). 




(4.11) 
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Consider the sum of the four diagrams for Wji shown in Figures 14, 15, and 16, paying 
particular attention to the r, a arrow flow. These diagrams correspond to those of Fig. 8, 
but the causality (or r,a) structure is clarified by the arrows. Indeed, examining these 
diagrams, we see that the soft fermionic line is either fully retarded or fully advanced, as 
indicated by uni-directional arrow flow along soft fermionic lines. As in the leading order 
case, this is a consequence of the fact that soft fermionic rr propagators are suppressed. 
The cut fermionic lines (those without arrows) are hard, and are only a function of p~ 
and not as is typical of an eikonal approximation. Thus, the p'^ integration is over a 
fully retarded or fully advanced function, and we are again free to deform the p"*" contour 
as in Fig. 13. After this deformation p^ is everywhere large (albeit complex) relative to 
P~ ,P±, Q, and we are free to expand the integrand at large p~^. 

The leading contribution in this expansion should be (p"*")^ (arising, for instance, from 
the soft propagator width, which first arises at NLO and should give rise to precisely 
such a {p'^)^ contribution); and the next order should be {p'^)~^- Higher orders are sup- 
pressed and need not be considered. The {p'^)^ term gives rise to a pure linear divergence 
/ cPp±d'^Q J dp'^ (p^)^ F {Q , p±) . A linearly divergent p~^ behavior in the p~^ ~ gT region 
corresponds to a leading-order behavior for p^ ~ T, so the linear divergence must ap- 
pear as a ^'''-independent small-p'*' limit of a hard-p''', leading-order contribution. There 
is precisely one such contribution, namely, the small p^ limiting behavior of the collinear 
region found in Eq. (3.11) of Subsec. 3.1.1. Therefore the 0{g), {p'^)^ behavior must be 
precisely Eq. (3.11), which was already included in the treatment of the collinear region 
and so should be subtracted to avoid double counting. 

Next we consider the subleading {p^)~^ behavior. In the last subsection we saw that 
such behavior arose already at the leading order, and that its physical interpretation was 
as an asymptotic thermal mass. While it is not obvious, it is at least not surprising that 
the subleading contribution should be precisely a shift to Eq. (4.9) in which is replaced 
by + 6m1^, as defined in Eq. (B.27). If we make this replacement and then expand to 
linear order in Sm^, we find 



p]_ + m'i, + 5ml^ p\ + ml^ °° \p\ + ml^ {p]_ + m^^f 

= + — . (4.12) 

p'i + ml^ °° {p\ + m2^)2 

The first term corresponds to the leading-order result and should be subtracted off; the 
second term is a true NLO correction. That is, we expect that the soft contribution at 
NLO should be 



(2^) 



3 d6T^ 



ST 2 2 , npik) f d p_L 2 P± 1o^ 



soft 



s 



k J (27r)2 - {pl + mir- 



Note that, like the leading-order term, this will also give rise to a logarithmic large-p_L 
divergence, which must be balanced by some logarithmic behavior at larger pj_. In this 
case the corresponding logarithmic behavior will be found in the semi-collinear region. 
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In conclusion, we expect three contributions from the soft region; the leading-order 
contribution which should be subtracted off, the infrared limit of the collinear contribution 
which should also be subtracted, and Eq. (4.13). The argument supporting this result is 
not rigorous, so we need to proceed with the actual evaluation of each diagram, making full 
use of the p'^ contour deformation technique. Since p"*" can always be taken as large, we will 
actually not need the vertex HTL's at all, and each diagram will become an expansion in p'^ 
as described above. We then sum the diagrams to get a gauge invariant total, whereupon 
we can perform the J d'^Q integral to find that we indeed get exactly the behavior described 
above. 

We believe that the very simple form of the large p"*" expanded result can be made rig- 
orous and understood physically in terms of eikonalized dipole propagation in the medium, 
and that this method can then be extended to other problems such as gluon radiation; we 
plan to return to this topic elsewhere. 

4.3 Soft diagrams 

We now turn to the diagram-by-diagram evaluation of the soft diagrams. The purpose of 
the next three subsections is to support with a concrete calculation the arguments of the 
last subsection, yielding in the end the same result presented there. We will concentrate on 
computing Wr; Wa is trivially related, corresponding to a contour in the other half-plane. 

The diagrams that contribute to the soft region at NLO are those in Fig. 8. We 
parametrize the different contributions of the soft NLO diagrams by 

diags. -1 

= ^ {Ws^ + + - W<^,,) (4.14) 

soft ^'^ 

where includes contributions arising from the two first diagrams of Fig. 8, i.e., from 
diagrams where the additional soft propagator gives a self-energy correction to the soft 
fermion. arises from the diagram where the hard fermion receives a self-energy cor- 

rection and gives the contribution of the remaining "cat eye" diagram. The last term 
^subtr arises from a mistreated kinematical region in the leading order calculation where 
the momentum in one of the hard loops in the resummed HTL propagator becomes soft, 
shown in Fig. 9. 

For each diagram there are several r, a assignments of the propagators and vertices that 
may contribute to the diagram. In practice though, most of these are suppressed by powers 
of g and give subleading contributions. The power-counting in the pure glue theory has been 
worked out in [28] . Here we extend the power counting to a theory with fermions. Because 
Gr oc l/p^ while Sr oc 1/p, the soft retarded gauge and fermionic propagators scale as 
l/{gT)'^ and 1/gT respectively. The symmetrized rr gauge propagator is proportional to 
ngijp) -|- 1/2 ~ T /p^ ~ 1/5 and is therefore enhanced by an extra factor of 1/ g; contrarily, 
the fermionic rr propagator is proportional to —np{p^) -|- 1/2 ~ p^ /^T ~ g and so is 
suppressed by an extra factor of g. Therefore in practice one always needs diagrams 
with the maximum number of soft rr gauge boson propagators, but no soft fermionic rr 
propagators. These rules are summarized in Table 1. 



d^k 
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Propagator 


Scaling 


Soft, retarded gluon 
Soft, rr gluon 
Soft, retarded fermion 
Soft, rr fermion 


l/T 



Table 1. Power counting for the soft propagators. 

All bare vertices have an odd number of a indices; those with one a index are the same 
as the zero-temperature ones, while those with three a indices carry a factor of 1/4. The 
three-a vertices remove Bose stimulation factors by reducing the number of available rr 
propagators and hence lead to subleading corrections in the current context, and we will 
not encounter them. 

Fig. 8 suggests that we would need the explicit forms of the HTL quark-gluon vertex 
and the two-quark, two-gluon vertex. However, while the HTL vertices are the same order 
as the bare vertices when all momenta are ~ gT, if a momentum p entering an HTL 
vertex becomes large p ^ gT, the vertex becomes suppressed, relative to the bare vertex, 
by {gT/p)'^. We will see that sum rules allow all diagrams (or combinations of them) 
to be evaluated in terms of their large-momentum behavior; therefore the total effect of 
HTL vertices will cancel exactly at NLO and we will not need their detailed form in the 
calculation. Therefore we defer their treatment to App. E. 

4.3.1 The soft-soft self-energy diagrams: Ws 



K + P 




P + Q 



Figure 14. The retarded diagram for the soft self-energy insertion. The arrows indicate r, a flow. 
Internal lines without arrows are understood to be rr propagators. The direction of fermion flow 
is clockwise and fermion momenta are always oriented along fermion flow. Wherever the arrow of 
r, a flow is parallel (antiparallel) to fermion flow this gives rise to a retarded (advanced) propagator. 

We start with {K) , the contribution to the retarded correlator from the self-energy 
insertion on the soft line. We observe that deviations from eikonality on the hard line are 
now NNLO at the largest, being suppressed by P/K ~ g, thus effectively reducing the 
hard line to an integral over x'^ of an eikonal Wilson line in the same direction multiplied 
by f^. and the appropriate statistical factor, i.e., {np{k) — l/2)fj^6{ivk • D). 
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As in the leading order case, we evaluate the retarded correlator and apply the KMS 
relation of Eq. (4.1) to obtain the backward Wightman correlator. Therefore, the left 
external line is of type r and the right one is of type a. There are several possible choices of 
r, a assignments for the internal lines, but given the power-counting rules in Table 1, only 
one represents an order- (7 correction to the leading-order result. It is the one where the 
gluon is rr, thus receiving a.l/g Bose enhancement, and the hard fermion is rr too. This 
r, a assignment and its corresponding momenta are shown in Fig. 14. 

Enforcing eikonality on the hard line, Srr{K + P) = —^^(1/2 — np{k))2'n5{2vk ■ P), the 
retarded amplitude reads^ 

W,^{K) = -2e'Y.qydng'C,J J " ^-(^)) 2v^5(^;fc • (Q) 

xTr [f^SAiPh^SAiP + Qh'^SAiP)] (4.15) 

where the factor of 2 accounts again for the possibility of having either of the two lines soft. 
Let us recall that in the soft approximation G^^^{Q) = T /q^ p^y{Q) and that the retarded 
HTL propagators in Coulomb gauge are given in App. A. 

Performing the p~ integration in Eq. (4.15) over the (5-function and defining for con- 
ciseness B{k) = e? q^e^diig^Cii{l/2 — np{k)) we have 



,SA{Ph''SA{P + Q)YSA{P) 



(4.16) 

This expression is a fully advanced function of p"*" . Using the analyticity arguments intro- 
duced in the previous section, we can again deform the contour away from the real axis 
in the lower half-plane without encountering poles or branch cuts. Let us call Ca the arc 
going from — ie to — at = ^+ ^ 5T, Im(p+) < 0. The integrand simplifies 
dramatically along this integration contour: the result of the trace and the propagators 
can be expanded for large \p'^\. As we argued, we need only the terms up to order l/p"^. 
We then have 

W^{K) = -2B{k) 



(4.17) 



where Gt is defined as Gij{Q) = {5ij — qiqj)GT{Q) (see Eq. (A. 4)), terms proportional to 
• q^, that average to zero in the azimuthal integration, have been omitted, and 

_ pl + ml^ _ (p^ + q^)2 + 




* We remark that there is an extra subtlety for fermions, since Sr{—P) — —Sa{P)- To assign the 
right prescription to propagators one should consistently assign fermion momenta parallel to fermion flow: 
wherever the arrow of r, aflow is parallel (antiparallel) to momentum/fermion flow this gives rise to a 
retarded (advanced) propagator. 
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We have furthermore used the fact that along Ca the components of S{P + Q) become 



+ Si(F + <3)^5^, (4.19) 

where in obtaining Eq. (4.17) we have expanded the "+" component for small (JE'p+q, since 
there is no other pole for g~ = on the opposite side of the complex plane, which would 
cause a pinch singularity. This will no longer be true when evaluating the cat eye diagram. 

4.3.2 Soft leading order subtraction: VFsubtr. 

Let us now turn to the subtracted counterterm of Figure 9. The entire advanced HTL 
self-energy, which is given by the simple one-loop self-energy graph taking one of the two 
bare propagators in the loop to be rr and the other to be advanced, with loop momentum 
hard and external momentum soft, results in the well known Shtl ~ gT. The self-energy 
we have inserted in Fig. 14 is however 0{g'^T) by construction. This implies that we have 
already implicitly subtracted all of the 0{gT) HTL self-energy and we have to only worry 
about 0{g) regions in the calculation of the HTL self-energy, where the approximations 
taken for its derivation fail. The only such region is the limit where the gluon becomes 
soft, which clearly overlaps with the phase space of the calculation we have just performed. 
A certain care is then needed in subtracting only this part of the HTL self-energy. To this 
end we take Eq. (4.16) and replace G''^(Q) with G(°)^^(Q) = T/q^p^^\Q), the soft limit 
of the bare gluon propagator {p^^\Q) is the bare spectral density). For what concerns the 
fermion propagator, we replace Sa{P + Q) with the bare one and, following with the HTL 
approximation, we keep only Q at the numerator, i.e., 

H/,,btr. ) - J (2^)3 J (2^)4 p2 + g2 + 2P • Q + ze{pO + qO) 

(4.20) 

The attentive reader might think that we could be neglecting other 0(g) regions, since in 
general, for a soft particle, the one-loop self-energy is equal to the HTL self-energy plus 
0{g) corrections, such as those that would arise from the inclusion of at the numerator. 
However the sum rule approach leads us to evaluate the self-energy at ^ gT and very 
close to the light-cone, where the full and HTL self-energies agree (as long as loop momenta 
remain hard, of course) and those 0{g) corrections vanish. Indeed, we have explicitly 
checked that the inclusion of ^ in Eq. (4.20) leads to vanishing extra contributions, at 
least up to order l/p~^. 

We can then evaluate Eq. (4.20) in analogy with the previous calculation, expanding 
along Ca- We obtain 

where we have used the fact that the bare longitudinal spectral density vanishes in Coulomb 
gauge. 
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K + P + Q 




p 



Figure 15. The retarded diagram for the hard self-energy insertion. Conventions are as in Fig. 14. 
The effective one-loop propagator for the hard line is understood to be rr. 



4.3.3 The hard-soft self-energy diagram: Wh 

In this subsection we compute the contribution arising from the hard-soft self-energy dia- 
gram of Figure 15 which we denote Wh- The power counting requires again that the soft 
line be advanced and the hard line, considered as a one-loop propagator, be rr, as shown 
in Fig. 15, together with the chosen momentum assignments.^ The KMS relation gives 
Srr{K + P) = (1/2 -np{k + p)){Sr{K + P)- Sa{K + P)), as in Eq. (2.3). The one-loop 
retarded and advanced propagators are obtained by inserting the retarded or advanced 
self-energy in the corresponding propagator. When plugged in our diagram, the advanced 
term gives rise to a fully advanced loop and thus vanishes, since all poles are on the same 
side of the p~ complex plane. The expression for W^{K) reads 

W, (K) - -mi) J ^ J + (4-22) 

where the Dirac and Lorentz structures of the hard line are again those of an eikonal Wilson 
line, i.e. v^v^f. Furthermore, the Wilson line gives rise to the retarded eikonal propagators, 
which are functions of p~ and q~ only. We can thus exploit their independence on p'^ to 
deform the contour along Ca and expand Tr [^^S'^(P)] along it, yielding 



W^{K) = +AiB{k) [ 



dp+ [dp- f d^p^ f d^Q G++(Q) 



271 J 27T J {2ttY J (27r)^Vk ■ P - ie)^Vk ■ (P + Q) - ie) 

+iV(4.23) 



2(p+)2(p- - 6Ep - ie) p^ 

where the {p~ — 6Ep — ie)~^ in round brackets on the second line has not been expanded 
for small 6Ep due to the presence of a pinch singularity with the double pole at p~ = —ie. 
We can then perform the p~ integration by closing the contour above, picking the residue 
from the pole in the first term in round brackets 

(A) - 2^B(k) J^^^ 2^ J ^2^^, J ^2^^, ^ ^ (4.24) 



^ The momentum assignments differ from tliose adopted in Fig. 1 and in the semi-collinear calculation. 
They amount to a shift, which will be undone in the end, when we will perform the subtraction of Eq. (3.11). 
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We now observe that the Q integration is free of pinch singularities when 5Ep goes to zero, 
so that we can safely expand the denominator. Furthermore the change Q ^ —Q (recall 
that GrriQ) is even in Q) brings the resulting expression to be identical to some of the 
terms in Eq. (4.17), i.e. 

W,{K)- 2^B{k)J^^ 2^ J ^^^^^J (2vr)4(p+)25i?2 ['+q--^e) ' ^""'^'^ 
4.3.4 The cat eye diagram: Wc 




Figure 16. The assignments needed for the O (g) soft correction to the cat eye diagram. Con- 
ventions are as in Fig. 14; in particular, internal lines without arrows are understood to be rr 
propagators. 

Next we consider the cat-eye diagram of Fig. 16. We label its retarded amplitude 
W^iK). The two assignments contributing are shown in Fig. 16. For the assignment on 
the right we have operated a shift, so that the momentum flowing in the hard rr propagator 
is always K + P.^^ The amplitude obtained by summing the two assignments in Fig. 16 is 
then 

xG';,{Q)v^Tt [^,Sa{P + Qh'SAiP)] , (4.26) 

where we have used the fact that Tr [f^SAiP + Q)YSa{P)] = Tr [f^SA{P)Y Sa{P + Q)] ■ 
At this point we integrate over the 5- function in dp~ , i.e., 

..,,/^/|;^(^.^) 

x5G-(Q)t;^Tr [^,5a(P + Q)7"5^(P)] |p-=o-(4-27) 

We observe again that the resulting function is fully advanced in p"*", allowing for a defor- 
mation along Ca- We furthermore notice that the retarded eikonal propagator {vk-Q — ie)~^ 
introduces a pinch singularity in the Q integration, since S^{P + Q) turns into the form 
given in Eq. (4.19). In order to make the pinch explicit we rewrite the terms in round 
brackets above as 

W-iK) - 2«,., I'^l^ (^^ + . Q)) 

xe;;;.{Q)v''Tv y,SA(p + q)YSa(p)] (4.28) 

This is analogous to what was done for the hard-soft self-energy. 
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Upon deforming the contour to Ca we obtain 



W^{K) = 2iB{k) 



Ca 



'2^ 



dq^(fq± 
(27r)3 



dq 



2tt {p+Y5Ep{q- - ie) 



pi + p±-q± ^++ 



\p+)HE^8E^+^ 
.g+('?i + P±-q±)G5TX'Z+,'Z±) 



p-^ 



p± • q± + Pi - ^Epp^ 



(4.29) 



where we have again dropped terms proportional to p_L-q± in the trace wherever they would 
have vanished in the angular integration. We furthermore observe that the terms propor- 
tional to q'^ /p^ on the second and third lines vanish upon integration, since GrriQ~^ , Q±) is 
an even function of q~^ . 

4.3.5 Summary and result 

We can now sum Eqs. (4.17), (4.25) and (4.29) and subtract the counterterm given by 
Eq. (4.21) to obtain the NLO retarded amplitude W^{K) = W,^ + W^/f + - W^^btr.- 
Most of the terms proportional to {q~ — ie)~'^ cancel, yielding 



W^{K) = -2iB{k) 



d'p± 
(2^)2 



Ca 



dp+ 
'2^ 



d\ 



{2'Kf 



pi 



{p+)HEp6Ep+^ 



(4.30) 




9i 



2p+{q- 



{p+)HEl 

where the last term is the subtracted counterterm. Furthermore, in dealing with the 
contribution from Eq. (4.17), we have dropped the term proportional to p±-q_L in ^E'p+q = 
5Ey, + {qi + 2p^ • q_L)/(2p''~). Exploiting the even nature of the gluon rr propagator we 
can rewrite the first term in round brackets on the second line in terms of 5{q~), yielding, 
after some rearrangements, 



W^{K) = -2iB{k) 
dq^ 



fp± 
(2^)2 



Ca 



2lT 



d\ 

(2^)3 



iG++(g+, g_L) ( p]_ + P± • q± 



{p+f 



5EI 



dEpdEpj^q^ 



+ 



pi 



27r 2{p+f5El 



{q~ - ie)2 



+ G?^(Q)|1 + ^ 




2G. 



(0) rr 
T 



iQ) 



(4.31) 



The first line is independent of p"*", since 5Ep, SEp^q l/p'^ ■ We already encountered the 
dqz integral in Eq. (3.4) (see also App. B.2): 



Tmi 



qi{qi + ml)' 



(4.32) 
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The second line is proportional to We show in Eq. (B.9) that the Q structure, 

without the subtraction of Gj, , can be identified as the two-dimensional condensate Zg 
defined in Eq. (B.l) when expressed in Coulomb gauge; 



(2vr)4 



Id 




1 




Zr, 



(4.33) 



( 0) TV 

The subtraction of Gj, precisely removes the leading-order contribution, leaving the 
NLO correction to Zg, ^Zg. [To see this, note that in Coulomb gauge the longitudinal 
spectral density vanishes and the transverse one is proportional to 5{Q'^).] In the UV, the 
bare and resummed propagator become identical, up to suppressed 0{m%/Q^) corrections, 
so that we can safely integrate up to infinity; therefore 6Zg is finite. We evaluate it in 
App. B.3, finding 6Zg = — Tm£)/(27r), a result originally due to Caron-Huot [22]. 
Plugging everything into Eq. (4.31) we obtain 



W^iK) = B{k) 



(2^)2 



'2^ 



(2^)2 (p 



+\2 



5EI 
Trrir. 



Trrir 



{p+f5El 27r 

We can now perform the straightforward dp'^ integral along Ca, leading to 



(4.34) 



W^{K) = B{k)T 



d^p± f d'^q± 



rrir 



(2^)2 y {27ry7r{pl+ml^)ql{ql+ml) 



pI 



pj + p± • q± 

(p± + q±)^ + 



-B{k)T 



d'^P± 



pi 



mj 



(27r)2(p2 +^2^)2 ^ ' 



(4.35) 



where we observe that ^'''-independent terms give a linear divergence. We can now plug 
this result in the KMS relation (4.1) to obtain the corresponding Wightman amplitude and 
from that the soft contribution to the rate, which reads 



diags. 



soft 



16/i+ f d^p^d'^qA 



(2vr) 



qi{qi+m%) 



P]_ + 



(p±+q±)^ + "^^ 



+ 



Y,^qie^dRnp{k) 
k 



p\ 



(2^)2(p^ + m2,)2' 



(4.36) 



where we have shifted the terms on the second line of Eq. (4.35) as in footnote A}^ We 
recall that the expression for the NLO correction Sm^ to the asymptotic mass is given by 



^^Had we used the momentum assignments of Fig. 1, we would have obtained the result directly in this 
form, see footnotes 9 and 10. 
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Eq. (B.27). This shows clearly how the contribution on the second line is nothing but what 
would have been obtained by substituting — )• + dm^ in the leading-order result 
(4.9) and expanding in g, as we anticipated. 

The first two lines in Eq. (4.36) are equal to Eq. (3.11), which is the soft-p^ limit of 
the LO collinear region we have analyzed in Sec. 3.1.1. There we concluded that Eq. (3.11) 
was to be subtracted from the rate obtained here in the soft region. Doing that removes 
the linear divergence and yields 



soft 



diags. 



soft 



(27r)3^^^ 



d^k 



subtr. 



soft 



k 

-^^(fc) 2 In 



Pi 



(27r)2 {p\ + m- 



NLOn2 



+ 1 



moo ) 



2 ^2 • 

OO/ 

(/zNLO)2 



NLO-)2 



+ mi 



(4.37) 



where the dp_i_ integration has been cut off at fj^^ . According to the analysis in Sec. 2.2 
as summarized in Fig. 10, we expect that the UV divergence will be removed by an IR one 
in the semi-collinear region, which we analyze next. Hence obeys g^T"^ <^ [fi^^)"^ <^ 

The result justifies the simplified approach in Subsection 4.2. At NLO, the contribution 
of the soft region is precisely the soft limit of the collinear contribution, the leading-order 
soft contribution modified by the shift 



''OO 



+ ^1^001 ^"^^ nothing else. 



5 The semi-collinear region 

As we have seen in previous sections, we must treat separately the region Q ~ gT and 
gT ^ pj_ <^ p'^ . Specifically we consider ~ gT^ while p+ ~ T, which means that the 
angle between p and k is small but not as small as in the collinear region; hence we will 
refer to this as the semi-collinear region^^^ i^j^ig case, one must compute the diagrams 
shown in Fig. 17, and then apply systematically the expansion k,p'^ ^ p± 2> q±,q^- 

Actually we have already evaluated these diagrams using the collinear expansion, since 
it is precisely these diagrams which give rise to the linear-in-collisions expressions we found 
in Subsec. 3.1.2. In particular, Eq. (3.13) was derived by making an expansion in ^ p±, 
and it still applies, with one proviso. In evaluating the collision strength C{q±) in Eq. (3.4), 
we treated p± ^ q_L ^ gT, leading to 6E ~ g^T. This let us neglect 5E when working out 
the kinematics of the gauge bosons. But if ~ gT^ then 5E ~ gT and this is no longer 
permissible. In particular, when we put (P + K) and {P + Q) on shell (see the cuts in 

^^With the momentum assignments in Fig. 1, a semi-collinear P and a soft Q imply that the momenta 
flowing through all fermionic lines {P + Q, K + P, K + P + Q) are semi-collinear. 
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Fig. 17), we find 
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where the 0{^) correction comes from p_L • q_L and always vanishes in the angular inte- 
grations. Therefore we must re-derive Eq. (3.4) with these somewhat different kinematics. 
A straightforward computation shows that the quantity 
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physically interpreted as the momentum diffusion coefficient and present in Eq. (3.13), 
should be replaced with its finite 6E generalization. 
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This expression reverts to Eq. (5.3) in the limit 5E — )• 0. Physically it represents the result 
of integrating over the cut gluon line in the scattering diagrams shown in Fig. 3, treating 
that line as soft. In the case 5E = kinematics force the gluon line to be in the Landau 
cut, but for 5E 7^ it can also be on-shell (on the plasmon pole); therefore both processes 
will contribute, so Eq. (5.4) will smoothly go over from the collinear splitting rate to the 
hard scattering rate as we increase p±_ and hence 5E. 

This (5£^-dependent momentum diffusion coefficient can be evaluated using Euclidean 
methods, see App. B.2, with the result, see Eq. (B.17), that 
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Unfortunately, performing the integral using Euclidean methods obscures what part arises 
from the Landau cut and what part arises from the plasmon pole. 
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Using Eq. (5.5) rather than Eq. (5.3) in Eq. (3.13), we find 
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The d^q± integration here is UV divergent, but the divergences will be removed when we 
subtract the two 0{g) regions of the leading order calculation that overlap with the semi- 
collinear phase space and which have already been included in those calculations. One such 
region was discussed in Sec. 3.1.2; when we evaluated the collinear splitting rate, we inte- 
grated over this phase space region but made the approximation that 5E — )■ 0. Therefore 
we should subtract q{6E = 0) from q{6E) used above. The other contribution we must sub- 
tract is the semi-collinear part of the phase space of the hard 2 2 contribution. The con- 
tribution we already included when we performed the hard 2^2 calculation corresponds 
to treating the gluon as free and on-shell. Therefore, we should also subtract from q{5E), 
its value obtained by using the free gluon propag ator, G^^(Q) o G(°)''''(Q) = T/q^p(^\Q), 
in Eq. (5.4). We will call this quantity g((5£')|bare- Explicitly, we find 
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where, as in our treatment of Eq. (4.33), we have used the fact that in Coulomb gauge the 
longitudinal spectral density vanishes and the transverse one is proportional to 6{Q'^) to 
simplify the integrand. We have used again the Euclidean techniques of App. B.l for the 
evaluation of the integral, although it can simply be evaluated in real time from the simple 
form of p^^\Q) = sgn(gO) 27r 6{Q^). 

Upon subtracting the two semi-collinear limits, i.e., 
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which is convergent in q^. The region where q± ^ gT represents a very small contribution 
to the integral, allowing us to extend the integration over the whole range. 

In order to evaluate the remaining integrals it is convenient to operate a kinematical 
distinction. Given our momentum assignments and the fact that p'^ ^ p^ + q^ , p'^ + k 
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p^+k, we can clearly see that the region where p'^{p'^+k) is positive corresponds, in terms of 
elementary processes, to having a semi-collinear quark or antiquark both in the initial and 
final state, i.e., to a Compton-like (for timelike Q) or bremsstrahlung-like (for spacelike Q) 
process. Conversely, a negative p+(p+ + h) is associated with the pair annihilation region, 
with a qq pair in the initial state. 

The integrand in Eq. (5.9) is symmetric under p'^ — >• —k — p"*", so the brem/Compton 
region is given by 2 dp'^ and the annihilation region by 2 /^^^g '^P 
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We have evaluated these contributions numerically using the same cutoff /i^^^ for the IR 
divergent p_L-integral as in the soft region (see Eq. (4.37)); the details are in Appendix. D. 
Summing up the contribution from brem/compton and pair processes, i.e., Eqs. (D.ll) 
and (D.21), we have the full contribution from the semi-collinear region. It reads 
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where the functions Cbrem/compton aiid Cpair are fitted by Eq. (D.6) and Eq. (D.17), re- 
spectively. Upon comparing this expression with the final result in the soft region, namely 
Eq. (4.37), we notice how the dependence on /i^^^ drops out of their sum. 

6 Results 

We can now collect all contributions and write the final result for the NLO calculation. 
Let us parametrize it as the sum of the leading-order result and its 0{g) correction 
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where the LO result is given by Eq. (2.8) and the 0{g) correction can be obtained by 
summing the collinear contribution, Eq. (3.22), the soft one, Eq. (4.37), and the semi- 
collinear one, Eq. (5.11), yielding 
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The dependence on the regulator /i^'"^ cancels in the sum of the semi-collinear and soft 
regions, as anticipated. The function Csoft+sc{k/T) is obtained by summing the non- 
logarithmic terms in the semi-collinear and soft contribution. It reads 
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Figure 18. Left: the function, C{k/T), parametrizing tlie photon emission rate for Nc = Nf = 3 
and as = 0.3 (see Eq. (6.8) and Eq. (2.9)). The full next-to-leading order function (Clo+nlo) is 
a sum of the leading-order result (Clo); a- collincar correction (JCcoii), and a soft+scmi-coUincar 
correction (i5Csoft-i-sc)- The dashed curve labeled Clo + (5^0011 shows the result when only the 
coUinear correction is included, with the analogous notation for the Clo + (^Csoft+sc curve. The 
difference between the dashed curves provides a uncertainty estimate for the NLO calculation. 
Right: the same as on the left but for larger k/T. 



Finally, we recall that A{k) and k are given in Eqs. (2.9) and (2.10) and (5m^/m^ = 
— 2?Tio/(7rT), as given by Eq. (3.25). The correction C^^n is intrinsically nonabelian, but 
dm^/m^ is nonvanishing in an Abelian theory. 
We now plot our results. Let us define 
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Given those definitions, it then follows that 
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= 0.3 and Nc = 
10, the collinear 



In Fig. 18, we start by plotting the function Clo+nlo(^/^) for cxg 
Nf = 3. In the phenomenologically interesting momentum range, k/T 
and semi-collinear-l-soft corrections largely cancel, leading to a small positive correction 
of order ~ 15% (Fig. 18(a)). At large momentum, k/T > 20, the LO and LO-hNLO 
curves cross and the NLO correction turns negative (Fig. 18(b)). We believe that the 
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Figure 19. The functions C{k/T) for A^c = 3, iV/ 3 as in Fig. 18, but for as = 0.05. 
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Figure 20. Left: the differential rate dT^/dk relative to the leading order rate as a function oi k/T 
(or cquivalently Clo+nlo/C'lo)- The full next to leading order rate (LO+NLO) is a sum of the 
leading order rate (LO), a collinear correction (coll), and a soft+semi-coUinear correction (soft+sc). 
The dashed curve labeled LO+coll shows the ratio of rates when only the collinear correction is 
included, with the analogous notation for the LO + soft+sc curve. The difference between the 
dashed curves provides a uncertainty estimate for the NLO calculation. Right: the same as on the 
left but for larger k/T. 



large cancellations we observe are rather accidental, and one should thus consider the 
curves Ci^o{k/T) + 5Ccon{k/T) and Ci,o{k/T) + 5Csoit+sc{k/T) as upper and lower limits 
respectively of an "uncertainty estimate" of the NLO calculation. 

In Fig. 19 we plot Clo+nlo(A;/T) and Ci^o{k/T) for = 0.05, and A^^ = 3, iV/ = 3. 
For the smaller coupling constant the NLO correction is always negative and rather flat, 
and the magnitude of the two largely canceling contributions is also significantly smaller 
than in the previous case. 



-38- 



(LO + soft+sc)/LO _i 



0.5 - 



0.5 - 



as=0.05 



as=0.05 











2 



4 



6 



8 



10 12 14 



10 20 30 40 50 60 70 80 90 



k/T 



k/T 



Figure 21. The differential rate, dTj/dk, relative to the leading order rate as described in Fig. 20, 
but for as = 0.05. 

In Figs. 20 and 21 we plot the differential photon emission rates dV^/dk relative to the 
leading order rate, (L0+ NLO)/LO, for two different values of the coupling constant. The 
reasonable, but somewhat ad hoc, "uncertainty estimate" described above can be inferred 
from the difference between the upper and lower dashed curves, which include either the 
collinear or the soft+semi-collinear correction, but not both. 

For the largest coupling, = 0.3, NLO corrections are modest and positive, although 
the "uncertainty band" is rather large - of order 50% (see Fig. 20). At intermediate 
coupling, Os = 0.15, the cancellation between the collinear and semi-collinear+soft contri- 
butions is quite dramatic, causing the LO+NLO result to be within a few percent of the 
LO rate (not shown). Nevertheless, the uncertainty band remains rather large - of order 
40%. Finally, at the smallest coupling Og = 0.05, the (LO+NLO)/LO ratio is somewhat 
larger than at intermediate coupling, but with a considerably smaller uncertainty band 



7 Conclusions 

We have computed the photon production rate to NLO of an equilibrated, weakly-coupled 
quark-gluon plasma. The contributions to the LO rate can be divided into distinct kine- 
matical regimes — the hard, soft and collinear regions. The contributions arising from the 
hard and the soft regions have logarithmic sensitivity to the details of how the kinematical 
regions are divided. However, this dependence cancels in the sum. At NLO the soft and 
collinear regions receive 0{g) corrections, and a new "semi-collinear" region starts to con- 
tribute here. We have dealt with the collinear region in Sec. 3, with the soft region in 4, 
and with the semi-collinear region in 5. 

The collinear regime is affected by the LPM interference of multiple scatterings through 
the integral equation (3.1). As we showed, computations are most easily performed in 
impact parameter space and the resulting 0{g) perturbation to the LO result is given 



(Fig. 21). 



in Eqs. (3.22). Furthermore two 0{g) regions of the leading-order result, the soft region 
and the semi-collinear region, are identified and the asymptotic behaviors are derived in 
Eqs. (3.11) and (3.13). The treatment of these 0{g) regions via the integral equation 
is incomplete, so we must recompute the contributions in these regions more carefully, 
subtracting off what has already been included in the collinear part of the calculation to 
avoid double counting. 

In the soft region, we relied on the fact that causality dictates analyticity properties 
for the n-point functions. Since the soft fields are lightlike separated as a result of the 
effective eikonalization of the hard fields, these analyticity properties lead to a tremendous 
simplification of the calculation: the soft bosonic correlators become the correlators of 
the 3D Euclidean theory, as pointed out in [21], whereas for fermions one can deform the 
integration contour of the non-vanishing light-cone momentum away from the real axis 
towards infinity, yielding a tractable expansion in inverse powers of that momentum. With 
these technical developments the soft contribution to the NLO rate is entirely analytic 
and given by Eq. (4.37), once the aforementioned subtraction of the collinear limit is 
performed. The result is incredibly simple: the NLO correction arising from the soft 
sector is the leading-order soft contribution modified by the shift — t- + (5m^, and 
nothing else. The contribution arising from the soft sector depends logarithmically on the 
UV regularization of the p^-integral; this dependence cancels against a corresponding IR 
divergence in the semi-collinear region. 

In the semi-collinear region, the contribution from the soft gauge fields factorizes into a 
light-like separated, two-point correlator given by Eq. (5.4). This correlator is related to a 
momentum dependent transverse momentum diffusion coefficient experienced by a particle 
whose momentum obeys the semi-collinear scaling. We use again Euclidean techniques to 
evaluate it, obtaining a compact analytic result. Performing the remaining integrations 
(the last of which only numerically) and upon subtracting the appropriate 0{g) regions of 
the leading-order result to avoid double countings, the contribution from the semi-collinear 
region is given by Eq. (5.11), which, as expected, shows an IR log divergence which combines 
with the soft contribution to make the final result finite and cutoff-independent. 

In Sec. (6) we combine all contributions together in the full NLO rate and plot it 
for different values of the coupling. The striking feature is that the NLO correction is 
composed of two largely canceling contributions. The positive one arises from the collinear 
region and the negative one from the soft and semi-collinear regions. For = 0.3, 3 colors 
and 3 light flavors, in the phenomenologically interesting momentum range, k < 15T, each 
contribution is a 50 — 75% correction, but their sum is just a 10 — 20% positive correction 
(see Fig. 18(a)). At higher momenta, the NLO correction turns negative, and the NLO+LO 
curve crosses the LO result (see Fig. 18(b)). For smaller values of the coupling, as = 0.05, 
the NLO correction is negative for all momenta (see Fig. 19). 

For these reasons we believe this cancellation to be largely accidental, and thus consider 
the two separate contributions as the upper and lower bound of an uncertainty estimate for 
the NLO calculation. In Figs. 20 and 21 we plot the (LO+NLO)/LO ratio of photo-emission 
rates together with the uncertainty band for the two values of the coupling. 

From the phenomenological point of view, the 0(20%) correction itself in the relevant 
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region for = 0.3 does not alter qualitatively the current analyses [30]. On the other hand 
our 0(100%) uncertainty band gives a first estimate on the reliability of the perturbative 
calculation. Going to NNLO one would also encounter UV vacuum divergences and the 
associated running coupling, whose scale setting introduces another possibly large error 
band. 

From the theoretical point of view, we believe that the main result of the present work 
lies in the developments related to the description of soft fields coupled to eikonalized hard 
fields. This progress opens new possibilities towards the calculation of other transport 
coefficients, such as the shear viscosity, at next-to-leading order. Furthermore, we believe 
the simple form of the NLO soft region can be understood more transparently in terms of an 
effective description of dipole propagation. Such a picture could also allow a factorization 
of the non-perturbative magnetic sector. We plan to return to these issues elsewhere. 

A clear extension of this work would be the NLO treatment of gluon radiation, following 
the generalization from photon to gluon radiation at leading order in [31]. Other possible 
extensions include the calculation of the NLO rate in = 4 SYM. The leading-order 
calculation at weak coupling was done in [32], together with the strong-coupling result. 
This could shed more light on the transition between the two regimes. 

Analogously one could apply the methodologies we have developed to similar calcu- 
lations for the thermal production of light-like particles, which could be of relevance for 
cosmology and whose rates are known only to leading order. Examples are ultrarelativistic 
right-handed neutrinos (see [29] for the LO rate), axions [33], saxions [34], axinos [35] and 
gravitinos [36]. 
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A Hard Thermal Loop propagators 

In this section we detail our conventions for the HTL propagators. Fermion propagators are 
most easily written in terms of components with positive and negative chirality-to-helicity 
ratio. The retarded fermion propagator reads 

Sn{P) = h;S:^{P) + h-S],{p), (A.l) 
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where 
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where the upper (lower) sign refers to the positive (negative) chirahty-to-hehcity com- 



ponent. The projectors are = (7° =p 7 • p)/2 and ml 
asymptotic mass squared. 

Gluons are described in the strict Coulomb gauge by 
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The other components of the propagators in the r, a basis can be obtained through Eq. (2.3). 



B Gauge invariant condensates 

During the calculation, we encounter several condensates that can be written as integrals 
of correlators separated by a spacelike or lightlike separation: 
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Here coordinates are written as triples {x~^,x , x±) with x = {t — z) and x"*" = (t + z)/2, 
so that t = z = x~^ when x~ = and X ■ P = x± • p± — x~p'^ — x^p~ . U{xi; X2) is the 
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Wilson line connecting the point xi to the point X2, in either the adjoint representation 
(Ua) or the representation of the fermions {Un), and W is the Wilson loop with a transverse 
segment of extent x± and a lightlike segment of extent x"*" in the {t,z) directions [37]. For 
Zg and Zf we have employed rotational invariance and chosen the light-like vector v to be 
Vk without loss of generality. 

The condensates Zg and Zf are related to the bosonic and fermionic hard thermal loops 
[38]. They describe how fluctuations in the gauge fields {Zg) and the fermionic fields (Zf) 
can influence the propagation of a fermion moving through the plasma at nearly the speed 
of light. Therefore they determine the dispersion correction of hard p'^ S> gT excitations; 
they are valid both at leading and at next-to- leading order [22]. And C{x_i_) and 
are related to the scattering processes in the medium; C{x±) arises when treating collinear 



splitting and is discussed in Sec. 3, while 
regime, in Eq. (5.4). They are related; 
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In this appendix we show how these condensates are most conveniently computed by 
using Euclidean methods developed in [21]. But first let us set up their calculation via real- 
time techniques, so we can see how expressions, encountered in the main text, do indeed 
correspond to these condensates. Except for C{x±), we only encounter the condensates at 
lowest order, where the Wilson line is set to unity and the field strengths take their abelian 
form. Working to this order, consider first Zg. Recall that 
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Applying this to Eq. (B.2) and performing the x integral, one finds (ignoring the difference 
between G'' and Grr, which is higher order in the soft region) 
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This is the same as the expression encountered in the "soft part" of the calculation, 
Eq. (4.33), except that there the leading-order behavior is to be subtracted. 

The calculation of q{5E) is analogous, except that the x'^ integral produces a delta 
function fixing q~: 
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This is identical to Eq. (5.4). 

It would be possible to compute these Minkowski-domain expressions explicitly using 
the sum rule approach; but we find it simpler and more instructive to compute them via 
Euclidean techniques. 
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B.l Relation to Euclidean functions: Simple derivation 

Caron-Huot has shown that n-point correlation functions, where all fields lie on a spacelike 
hypersurface, can be carried out by Euclidean techniques [21]. The null correlators we need 
can also be computed provided that they are free of collinear singularities, which they are. 
The proof presented in Ref. [21] is rather complex and technical. Here we will present a 
much simpler derivation, also due to Caron-Huot^^, which works for two point functions. 
In practice this is all we need, except for the NLO evaluation of C(xj_). 

Consider the ordering-averaged correlator of some operator (such as the field strength), 
Grr{x^,^) with \x^\ > \x^\. (Since the separation is spacelike, operators commute, and 
therefore Grr equals G^, Gp, or G^.) Write it in terms of its Fourier representation 

Grrix',^) = Jdc^l dp,d'p^e'(^'P'+^^-^^-'^'^'^ Grr{co,p,,p±), (B.ll) 

and use 

Grr{oJ,p) = (^Ugiuj) + ^ {Gr{u!,p) - Ga{u},p)) = (^Ugiuj) + ^ p{(^,p), (B.12) 

and define p^ = p^ — {t/x^)oj: 

G,,(a;°,x) = jdujj #,dVe'("'^"+"^ P^^ (ns{oj) + ^^ p{oj,p + u;ix'' /x'),p^) . (B.13) 

Now we perform the uj integration by contour methods. The retarded function in p{u,p) 
is related to the Euclidean function via Gr{ijJ,p) = —iGE{i^^,p), that is, by analytic 
continuation. This continuation is guaranteed not to encounter singularities in the Green 
function so long as the imaginary part of the 4-momentum remains timelike - since then, in 
some frame, the continuation is purely of the frequency. Since < 1, the continuation 

of Gr{uj,p — oo{x^ /x^),p±) in to will not encounter any singularity in the upper complex 
Lo plane. The advanced function similarly will be free of singularities in the lower plane. 
Therefore the only singularities encountered in continuing the frequency integration are 
those in the statistical function {ng{uj) + 1/2), which has poles at w = 2'KinT with n = 
(. . . — 1,0,1,...) and residue equal to T. Closing the contour around these poles, and 
renaming p^ to p^, we find"'^'^ 

Grr{x\x',X^)=TY, I ■^e'^^'P'+^^-P^^GE{i^n,Pz + iMx^/xn,P±) , Un = 2^nT . 

(B.14) 

In any case where we need to compute the soft gT contribution to such a correlator, 
one may drop the nonzero Matsubara frequency contributions; that is, we keep only the 
n = term in the sum. For this term, GE{i^n,Pz + i^nix^ /x^),p±) = Ge{0,PziP±) is the 
Euclidean correlation function of the 3-dimensional dimensionally reduced (Electric QCD 
or EQCD) theory. [EQCD is the 3D theory consisting of the spatial gauge fields Ai and an 
adjoint scalar descended from the temporal component of the gauge field, ^> = iAq. The 
scalar is massive, m| = m^; see for instance [39].] 



^^S. Caron-Huot, oral presentation at the Institute for Nuclear Theory (Seattle), 29 March 2012 

The pole at n = is an artifact of the separation of p, which vanishes for uj = 0, into Gr and Ga- 
The individual poles there can then be dealt with in a principal value prescription, for instance. 



-44- 



B.2 Application to Scattering 

First we apply this method to determine C{qj_) and q{5E) at leading order. There, 

CloC^x) = g^C^j dx+ j d2xxe-^^^-"^G++(x+,0,xx) , (B.15) 

where this position space expression makes clear that we are dealing with a spacelike 
separation and hence Eq. (B.14) is applicable, yielding 

= q^CbT ( 7^ —] + hard (n / 0) contributions. (B.16) 

\ql Ql+miJ 

The l/q'j_ and l/{q'j_ + m'^) terms are the contributions from the A^A^ correlator and 
the Ai^Aq correlator respectively. The nonzero Matsubara frequency contributions are 
suppressed, for ~ pT, by a power of and may be neglected. This result was first found 
by Aurenche Gelis and Zaraket [26], by a rather more complicated sum rule procedure, 
which involved canceling l/((?^ + m^/3) poles. The current procedure, originally due to 
Caron-Huot [21], avoids this complication. 

With this Euclidean framework the NLO (1-loop) corrections also become tractable as 
a computation within dimensionally-reduced EQCD [21]. We will return to this result in 
the next appendix. 

Next consider q{6E), Eq. (B.4). We work to lowest order, replacing the Wilson line with 
the identity and keeping only the two-point correlator of the ^-fields in the field strengths. 
We may again apply Eq. (B.14) and to find the infrared contribution we may keep only 
the n = term. This corresponds to the replacements Grr Ge, f dq^ /{2-k) — t- T and 
q^ —7- lirinT 0, q^ ^ q^- Writing the gauge field correlator in terms of its momentum 
space representation, we then obtain 



(2^)= 



d2 



(2vr)^ 



++ 



{0,q,,q±) + G^{0,q,,q^) 



x2Tr5{q^ + 6E) 



mlq\ 



+ 



1+- 



25E^ 



6E^ - 2q,5E 



{ql+6E^){ql + 5E^ + ml) ql+SE^' 



(B.17) 



B.3 Application to 5Zg 

Let us now apply the Euclidean formalism to compute Zg at NLO. Our starting point will 
be the position space expression of Eq. (B.2). Since the positions involved are lightlike 
separated, we may apply Eq. (B.14). 

The leading order contribution to Zg arises from hard p^ ^ T excitations; it therefore 
arises from a range of n values and is not easily established by Euclidean methods. However 
the leading order value is easy to determine by conventional methods, and equals 



Z. 



LO 



d^p , , r2 



(B.18) 
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For n 7^ the correlation functions receive 0{g'^) loop corrections, so any corrections 
from n 7^ are 0{g'^) and therefore beyond our current precision goal. To find the 0{g) 
NLO corrections, we need to compute the contribution from n = modes and subtract 
the leading-order, unresummed-theory value. A straightforward evaluation of Eq. (B.2), 
replacing x'^ with Xz and using 3-dimensional Euclidean correlation functions, gives 



5Zg=T 



(2vr) 



-2 



, . ,„ + Grp{Q,qz,q±) [?> 1 - same, free propagators 

(g^ + ley 



(B.19) 

where the ie prescription in l/{qz + ie)^ arises from the boundary conditions on the z 
integration, zdz. Now the bare and interacting Ai correlators are the same at leading 
order, but the temporal mode Aq develops a (Debye) mass; (^0^0(9) = i whereas 
the bare value is —1/q'^. (The minus sign is because there is an i difference between the 
Minkowski and Euclidean A^ field, so A~^ = A^ + iA^.) Therefore only the part of the 
first term is not canceled by the relevant free version. Its subtracted contribution is then 

" ' y (27r)3 (g, + iey W+ml^ q^J- 27r ' ^^'^''^ 

The correction is 0{g) relative to Eq. (B.18) and negative, representing a reduction in the 
thermal mass due to the screening of infrared gauge modes. 

No such 0{g) correction arises for Zf because the correlator directly involves fermionic 
fields which do not have zero modes. To 0{g) accuracy, 

B.4 Thermal mass at NLO 

We now apply these results for Zg and Zf to the fermionic effective mass. The thermal 
dispersion relation for a particle approaches Pq = p'^ + for p^ p gT, where 
is the asymptotic mass. The asymptotic mass is given by the real part of the fermion 
self-energy. In more detail we have [40] 



m 



2 -2pReS+(pO=p), S+(P) ^ ^Tr . (B.22) 



Q 




Figure 22. The diagram contributing to the asymptotic mass at leading and next-to-leading order. 

To the order of interest, one diagram contributes to the retarded self-energy, see Fig. 22. 
There are two r, a assignments; either the fermion is retarded and the gluon is rr, or the 
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gluon is retarded and the fermion is rr. The retarded hne, which in the NLO calculation 
always carries a large momentum, enforces the eikonality; the operator correlation in Zg 
(Zf) corresponds to the cut gauge boson (fermion) line. One easily finds that [22] 

ml = g^C^iZg + Zf) . (B.23) 

We just found that, at NLO, the condensates read [22] 

6 ' ^ 2tt 
12 

which, together with the well known result for the leading order Debye mass 



Zg = Z^^ + 6Zg, Z'f = t-^ sZg = -±^, (B.24) 
Zf = Zf + 5Zf, Zf = ^, 6Zf = 0, (B.25) 



3 

yields the fermionic at NLO in QCD,^^ 



l = ^{C^ + TRNf)T\ (B.26) 



rriocNLO = + ^^oo = 9 Cn \- —]j j ■ (B.27) 

C NLO collision kernel 

In the previous appendix we presented the definition of the differential collision rate C{x±) 
and found its leading order expression in transverse momentum space, C{q±). Its NLO 
expression in q± space has also been found [22]; writing Cnlo = C + 6C, explicitly 



mo+2^i^arctan^ - arctan 3^ 



arctan arctan 7 



,4 



+ 



27rq±{ql+ml) 2^g3 32^3 

rriD / 3 2 1 



(C.l) 



In order to use the collision kernel in Eq. (3.16), we need to Fourier transform this expression 
to find C'(b), preferably analytically but at minimum via accurate numerical integration. 

The integral at leading order can be performed by taking partial fractions and consid- 
ering e"^^'^-^ / {q'j_ + m%): 



d^l^b,, r dq2 1 ^ r dqi^k,, 1_ 



00 27r 2tt ql+ql+ m% J 2tt 2^ q{ + m\ 



1 dx _ Ko^bniD) 



2vr y^x2 - ml 27r 



(C.2) 



In [40] the asymptotic mass was computed numerically at NLO in the large- A'^/ non-abelian theory. 
The authors found a strong momentum dependence of their NLO result, in sharp contrast with Eq. (B.27), 
which however reproduces their average shift. It is our understanding that, in the context of a strict g 
expansion, the result of Caron-Huot [22] is correct for fermion momenta much larger than gT. 



-47- 



where, after performing the trivial q2 integration, we changed contours to pick up the 
discontinuity along the cut of the y^g^ + m|, function; x = Imq. Ultimately we need the 
same expression with mj, — t- 0, and with 6—7-0 for each (finite and zero rrio) case; 

r2 



q^TCb. lim ^ {Ko{ebemij) - Ko{bemo) - Ko{ebmo) + Koibrrio)) 



Cjoib) 1 



(Koibmi,) + 7E + ln(6mo/2)) . (C.3) 



We follow the same strategy for the terms in the NLO correction (C.l). We first 
perform the q2 integration (the direction orthogonal to b). This can generally be done by 
deforming the contour to pick up all poles and cuts. In every case the resulting integral can 
be done analytically. Then we perform the qi integration. In some cases this can be done 
analytically, in other cases the b = case is analytic but the finite b case involves an integral 
along a cut. Finally there are cases where the difference, 1 — e~^^ must be integrated along 
a cut. Some of these integrals remain numerical, but all converge exponentially and are 
small for large b. 

Without going into detail, the result (writing the terms in the same order as they 
appear in Eq. (C.l)) is 

mo 6C'{b) 



1 



2 A V J {z^- 1)3/2 



"(1 



967r2 Jo z \ z^ 



32 87r2 7i Vi^^ 87r2yo z 

lbm,o 1 



647r ' 87r2 



Ko{2bnia) - 2Ko{bmo) + In ^ - 7e + bmoKi{bmi,) - l] . (C.4) 

rriob J 

We have been unable to perform these integrals analytically. However, they are quite 
straightforward numerically, and they also have a simple behavior in the large 6m £, limit: 

mnSC'ib) 7bmn -Ovr^ - 58 + 72 ln(2) - 3(ln(6mo) + 7e) , ^ 

+ ^^^32 + 0(exp(-6m^)). 



g'^T^CnCA 647r 2887r2 

(C.5) 

This has the interpretation as the total scattering rate, IR regulated at a scale g ~ 1/6. 
While the leading order total scattering rate has a logarithmic divergence, the NLO has a 
linear divergence. Note that the corrections to this linear + constant + log behavior are 
exponentially suppressed. 

The small bmo behavior is 

ITT'D SC'jb) ^ -bmo qnlob'^ml 3 „ _ 37r2+10-41n(2) „s 
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Here ^nio is the NLO correction to the momentum broadening rate. To see this, expand 
Eq. (3.17) in small b and angle average, exp(ib • q^) ~ 1 — (b • q_L)^/2 ~ 1 — iP'q^ /A. The 
integral J cPq q^C{q) is what we usually mean by q; hence ^nio is the numerical coefficient 
on the NLO contribution to q. We have thus verified that the expansion performed on 
the expression in Eq. (C.4) reproduces the expected linear-in-6 and q behavior found by 
Caron-Huot [21]. 

Note that the small h behavior is actually negative. Our understanding is that this 
is actually correct. Basically, small b corresponds to large momentum transfers, a limit 
where our treatment of exchange momenta as gT <^ T breaks down. In [21], Caron-Huot 
showed how an opposite term arises in the soft limit of the hard contribution to q, leading 
to a cancellation. Had we used dimensional regularization to perform the integration of 
6^(/^Cnlo(9)/4, the linear term would simply have vanished. 



D Evaluation of the semi-collinear integrations 

In this Appendix we evaluate numerically the integrals appearing in Eq. (5.9). In order 
to match properly with the UV divergence in the soft region, which was regulated by 
a cutoff /i^^^ {gT ^ fi^'^ ^ \/ff^) the transverse momentum p±, we regulate the 
IR region with the same transverse cutoff. As we mentioned, we split the calculation 
into a bremsstrahlung/Compton contribution (p+ > 0) and an annihilation contribution 
{—k/2 < < 0), where we used the symmetry under p+ — t- — p"*" — k. We start with the 
former. 

D.l The bremsstrahlung/Compton contribution: r^^J^^^^)"™''* 

As we shall show, it is technically convenient to first introduce an intermediate regulator 
/i"*" for the dp^ integration, with gT ^ /i+ <C T, hence fi'^ » fi^^. In practice we divide 
the phase space in two regions, i.e. 

1. First we take p~^ > fJ-"^ , p± > 

2. We then consider the slice < < fi'^, p± > fJ^^ , where only the IR asymptotic 
behavior in p'^ needs to be considered. 

The dependence on cancels in the sum of the two regions. 

The p± and q± integrations on the second line of Eq. (5.9) yield 

1 f°° dp^p^ f d^qi_ \ mlql ml 



l-K y^NLo (5^2 y (27r) 



(g5_ + 5^2)(gi + '5^^ + m-2) g^ + m2 



= -Slp^(fc + P+)l + OM"°)^), (D.l) 

where we have used the fact that p"*" ^ /i^^^. The first term on the second line is the 
result one would obtain with vanishing cutoff. As we mentioned in the previous section, 
matching regions for momenta of the order of the cutoffs are equally described by the 
regions on either side of it and the dependence on the cutoff has to vanish at all orders. 
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Power-law terms can then be neglected and we can just plug the first term in Eq. (5.9), 
yielding 



(2vr) 



(1) 

brcm / compton 



A{k) 



dp 



np{k) 



kp+ {p+ + k) 



(D.2) 

This expression is IR log divergent. The logarithm can be extracted by adding and sub- 
tracting 9(T — p^)/(2p^) under the integral sign, so that 

rT 



ttT 



A{k) / dp 



1 



rrir 



2p+ 2ttT 



Aik)ln' 



and 



d^k 



(1) 



brcm / compton 

Cbrcm/compton(^/^) is defined as 



brem/compton I rp 



k 



mi 
'2^T 



dp^ 



A{k) 



I ) ~l~ C'brcm/compton 



(D.3) 



(D.4) 



np{k + p+){l - np{p+)) {k+p'^f + {p+f 



np{k) 



]^P+ _|_ /j;) 

e{T-p+) 
2p^ 



, (D.5) 



where we stretch the integral to since it is now finite, the difference being given by negligi- 
ble positive powers of fJ.'^ /T. For further convenience, we parametrize Cbrcm/compton(^/^) 
with an accuracy of 2% or better as 

t2 4 



Ch — 



-'brem / compton l-^ 

with 



V <i.e— + ( -0.12563 + ^ + 



n=l 



, (D.6) 



[di... di] = [2.29467534576455, -3.0183977101591, 1.2374580732449,0.486264291149683]. 

(D.7) 

Let us now turn to region 2. By imposing the cutoffs and expanding for p^ -^k we have 
(2) A rt^+ roo 



brcm / compton 



4 1'^ 

-A{k) / dp+ / 
T Jo Jf,] 

d'q^ 



P± 



(27r) 2{p+ySE^ 



(2vr)2 



mlql 



rrir 



{ql+5E^){ql+5E^+ml) ql+ml 



(D. 



where now 5E ~ p^/(2p^). By changing integration variable from dp'^ to d5E we have 



(2) 



"-A{k) 



dp± 



brem/compton ^ ^^^O (2^) ^i/(2/.+) P±SE^ 



d6E- 



d^ 



q± 



(2vr)2 



mlql 



rrir 



{ql+6E'^){ql+6E'^ +ml) ql+ml 



(D.9) 
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which is easier to integrate. Upon expanding the result in /i^ ^ fj^'^ we have 

brcm/compton 



m 



° A(k) hi ^"^""^ I O 



(D.IO) 



Let us remark that the 0{{fi^^)'^A{k))/{fj.~^T)) term obtained here cancels exactly with 
the one that would be obtained by keeping the suppressed ©((/i^^*-*)^) term in Eq. (D.l), 
confirming the correctness of our two-region matching procedure and the aforementioned 
cancellation of power-law terms in the cutoffs. 

Hence, the contribution from the brem/Compton region is 



brcm./Compt. 
semi— coll 



rrir 



2-kT 



A{k) 



In 



2Tm 



D 



(^NLO)2 



brcm / compton 



(D.ll) 



D.2 The pair annihilation contribution: rgcj^i„c,oii 



We employ the same two-region strategy for this process too. In principle we would have 
IR divergences in p+ at both endpoints and —k. However, using the symmetry at —k/2 
we discussed before, we can restrict the integration to (— /c/2,0) and worry about a single 
IR divergence. In practice, in region 1, plugging Eq. (D.2) in Eq. (5.9) and keeping only 
the /z^^^-independent term, we have 



(2vr^3^^^ 
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2^T 



A{k) 



-"^ , + n^(fc +p+)(l - n^(p+)) {k+p+f + jp+f 
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Upon changing the sign of the integration variable and restricting the integration to (0, k/2) 
we have 



(1) 



pair 



' ttT 



A{k) 



^^+Mk - P+) nAp+) [k - p+f + {p+f 
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(D.13) 



The IR logarithm can be extracted by adding and subtracting l/(2p~^) under the integral 
sign, i.e., 



27rT' 2/i+ 

With a slight rearrangement we can write the result from region 1 as 
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2p^ 



We can fit Cpair with a very good accuracy as 



CpUk/T) = -IE + Ei(-A:/2T) + I^^k/T), 



(D.14) 



(D.15) 



(D.16) 



(D.17) 
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with the function 



Ipair(rE) ~ 2ln{g{x)) + b, + + + ( a„xM e-^ (D.18) 

g{x) g{xy J 

where g{x) = e~^/^ + x and 

61 = -1.29715 , 62 = 1.38486 , 63 = -1.58046 , (D.19) 

and 



ao = 1.2014 , ai = -0.303456 , a2 = 0.00446236 , 03 = -0.0451118 . 



(D.20) 



Region 2, upon exploiting again the symmetry with respect to k/2, turns out to be identical 
to its brem/compton counterpart, yielding again a log that removes the dependence. 
The final contribution from the pair processes is thus 



(27r)3'^^^ 
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E The contribution from HTL vertices in the soft region 

In this section we analyze the contribution from the HTL vertices within the framework 
of light-cone fermionic sum rules we have introduced before. The analysis performed in 
Section 4.2 relied heavily on analyticity arguments; relations such as the KMS relation were 
employed to rewrite propagators or amplitudes in terms of fully retarded and fully advanced 
functions. This is particularly advantageous in the current analysis of the contribution 
of the HTL vertices, since, as we shall show, only the fully retarded/ advanced vertices 
are needed, i.e. only those with one external a line, which correspond to the analytic 
continuation of the Euclidean Hard Thermal Loops. 

As we mentioned, a full treatment of the HTL effective theory within the context of 
the r, a basis was carried out in [28] in the gauge sector only. There it was observed that 
hard loops with two external a lines are enhanced by a Bose factor of T/p^ ~ 1/(7, P 
being the external momentum, with respect to the ones with one external a line, which 
are instead the fully retarded functions obtained by analytic continuation of the Euclidean 
amplitudes. However, when attaching propagators to these loops, the enhanced ones can 
only be connected by ra propagators, which scale like l/{g'^T'^) (see Table 1), whereas a rr 
propagator can be attached to the standard, fully retarded ones. The rr propagator has a 
relative T/p^ ~ I/5 enhancement, thereby making the connected amplitudes of the same 
size. 

An opposite behavior is observed when the analysis of [28] is generalized to include 
fermions, as we have done. Consider for simplicity the HTL self-energy: the ra and ar 
amplitudes, which are the fully retarded and fully advanced amplitudes, corresponding to 
analytic continuation of Euclidean loops, scale like gT. The aa amplitude scales instead 
like g'^T, the suppression being due to Pauli-blo eking; this can easily be seen by noting 
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that the aa amphtude is hnked by the KMS relation to the difference of the ra and ar 
self-energies times (1/2 — npijp)) ^ /{AT) ~ g. 

Going to the effective quark-gluon vertex, one obtains that the raa amplitudes are all 
Pauli-blocked by a factor of g with respect to the fully retarded rra, irrespective of which 
particle is assigned the r label, whereas the aaa vertex scales like the rra, in agreement 
with the results of [41]. This would then give a more complicated set of power-counting 
rules than in the pure gauge theory. For the problem at hand, however, a limited number 
of hard loops is required: the requirements that the gluon be rr, in order to gain & 1/g 
enhancement, and that the soft quarks be ra or ar, in order not to be Pauli-blocked, imply 
that the only needed HTL vertices are of the simple rra and rrra type. For what concerns 
the cat eye diagram one just needs to replace the bare vertex on the soft quark line in 
Fig. 16 with the HTL vertex having the same r, a assignments. Regarding the soft-soft 
self-energy and tadpole diagrams, the three possibilities are depicted in Fig. 23. They are: 




Figure 23. The relevant diagrams and assignments for the soft, retarded self-energy with HTL 
vertices. The gluon is always rr. The notation for the HTL vertices follows the one in [28]: the 
double line is the eikonalized hard mode flowing in the loop. 



1. A single vertex connecting the gluon to the soft line is replaced with its fully retarded 
(or fully advanced) HTL counterpart, whereas the other vertex is kept in its bare form. 
One of the two possibilities is shown on the left in Fig. 23. 

2. Both vertices are replaced by HTL vertices, as in the middle diagram of Fig. 23. 

3. Only the two-quark, two-gluon rrra HTL vertex is relevant for the tadpole, as shown 
on the right in Fig. 23 

A detailed inspection shows that no other assignment can contribute at the same order in 
9- 

At the practical level, the quark-gluon effective vertex is obtained with this simple 
replacement 



47r {vi- P + ie){vi- {P + Q) + ie) 

in Eqs. (4.15) and (4.26). vi = (1, v;) is a lightlike four-vector, corresponding to the hard 
momentum in the loop (hence the label I), whose direction we integrate on. 
The contribution of the tadpole to the retarded, soft self-energy reads instead 



{vi-P + ief \vi-{P-Q) + ie 
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This expression contains a symmetry factor of 1/2. 

As we already remarked, the analytic structure in p"*" of the cat eye and self-energy 
insertion diagrams, which was crucial for the derivation of the sum rule, remains unchanged 
with the addition of the HTL vertices, so that we can still deform the integration contours 
at \p~^\ S> gT. In order to prove that the HTL vertices are irrelevant in the sum rule 
context, we then need to show that the amplitudes that include them go to zero on the 
arcs faster than l/p"*". Although the HTL vertices introduce two more powers of at 
the denominator, the fact that P is a spacelike vector and that P -\- Q might also be 
spacelike causes possible enhancements in the d^li angular integration. We thus proceed 
by performing the traces obtained by inserting the HTL vertices in the amplitudes Wc 
and Ws through the rules we have just introduced. One then obtains a basis of angular 
integrals of the kind 



where f{vi) is a set of functions obtained by contracting Vi with all other 4- vectors available, 
i.e., 



and products thereof. U = (1,0,0,0) is the four-velocity of the plasma in its rest frame 
and vi' is the hard loop velocity of the second HTL vertex, which arises in the evaluation 
of the middle diagram in Fig. 23. The contribution of the tadpole can also be related to 
this basis by differentiation with respect to and q^. 

The angular integrations are known in the literature and can be read from [42, 43]. Upon 
inserting the results back in the amplitudes Ws and Wc and expanding for large p~^ one 
obtains that all contributions containing one or two HTL vertices behave on the arcs as 
l/{p~^)'^ or smaller, and therefore do not contribute. 
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